EMBEDDED CONSTANT MEAN CURVATURE SURFACES IN EUCLIDEAN 

THREE-SPACE 



CHRISTINE B REINER AND NIKOLAOS KAPOULEAS 

Abstract. In this paper we refine the construction and related estimates for complete Constant 
Mean Curvature surfaces in Euclidean three-space developed in [10] by adopting the more precise 
and powerful version of the methodology which was developed in [14]. As a consequence we remove 
the severe restrictions in establishing embeddedness for complete Constant Mean Curvature surfaces 
in [10] and we produce a very large class of new embedded examples of finite topology. 



1. Introduction 

Critical points to the area functional subject to an enclosed volume constraint have constant 
mean curvature H = c ^ 0. We let "CMC surface" denote a complete, constant mean curvature 
immersed smooth surface in R 3 . The only classically known CMC examples of finite topological 
type were the round spheres and cylinders and more generally the rotationally invariant surfaces 
discovered by Delaunay in 1841 [3]. Hopf [8] proved that the only closed CMC surfaces of genus zero 
are the round spheres. In 1986 Wente constructed genus one closed immersed examples [28]. Using 
a general gluing methodology developed in [10,27], and using the Delaunay surfaces as building 
blocks, most of the possible finite topological types were realized as immersed CMC surfaces for the 
first time [10, 16]. However, no genus two closed examples could be produced with the constructions 
in [10, 16]. In [14] a systematic and detailed refinement of the original gluing methodology made it 
possible to construct genus two closed examples by using the Wente tori as building blocks. Since 
then, many other gluing problems have been successfully resolved by using this refined approach 
[6,7,12-14,16,18,19,29]. 

We discuss now the case of embedded, or more generally Alexandrov embedded, CMC surfaces. 
Let M 9: k denote the moduli space of all Alexandrov embedded CMC surfaces with H = 1, finite 
genus g, and k ends where two surfaces are considered equivalent if they differ by a rigid motion 
of M 3 . Using a reflection technique, Alexandrov [1] proved the only embedded £ € Mg,o is the 
round sphere. Meeks [26] proved the space A4 3i i is empty and that every end of £ G M. g ^ is 
cylindrically bounded. Motivated by [10,26], Korevaar, Kusner, and Solomon [21] showed that 
each end converges exponentially fast to a Delaunay surface and any £ £ A4 3i 2 is necessarily 
a Delaunay embedding. Kusner, Mazzeo, and Pollack [23] proved that Mg : k is a real-analytic 
variety. Moreover, for a non-degenerate £ 6 M g ,k there exists a neighborhood of £ in M. g ^ that 
is a real-analytic manifold of dimension < 3k — 6. Here strict inequality occurs when one quotients 
by the finite isotropy subgroup of £ in the group of Euclidean motions. 

Some further constructions of CMC surfaces have been carried out. Grofi-Brauckmann [5] used 
a conjugate surface construction to construct the surfaces in A4o,fc with maximal (fc-fold dihedral) 
symmetry. This includes those which have large neck size (the examples in [10] all have small 
neck size) . Mazzeo and Pacard [24] extended the gluing construction to produce CMC surfaces by 
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attaching Delaunay ends to the ends of a non-degenerate Alexandrov embedded minimal surface 
with finite total curvature, genus g, and k catenoidal ends. Mazzeo, Pacard, and Pollack [25] proved 
that a Delaunay end can be attached to a non-degenerate S G -M g ,k to produce a non-degenerate 
£' G M fl ,fc+i. 

In this paper we return to the general construction in [10] and refine the construction and esti- 
mates by applying the improved methodology developed in [14]. In particular we can now produce 
a large class of embedded examples of finite topological type. In [10] too much "bending" of the 
catenoidal regions was required to ensure a successful construction and this destroyed embeddness 
(but not Alexandrov embeddeness). It was possible to avoid the need for this bending only in a 
few cases of very high symmetry, more precisely only when a few Delaunay ends were attached to 
a central sphere and the symmetries of a Platonic solid were imposed. In the current construction 
no such bending is needed; therefore we can ensure embeddedness for all reasonable candidates. 
The examples we construct in Ai g ,k have 3A; — 6 continuous parameters as expected. Because no 
symmetries need to be imposed in our construction, we hope that it will serve as a model for ap- 
plying the general methodology in general settings without symmetries. The current construction 
can also be extended to higher dimensions [2]. 

We now provide a moderately detailed outline of the construction in the subsections that follow. 

The graph and the parameters. We begin the construction by designating a finite graph F 
that satisfies a few necessary geometric conditions and a flexibility condition. T consists of a 
finite collection of vertices, edges, and rays and to each edge or ray we assign a parameter f(e). 
The geometric conditions arise as a result of certain properties of CMC surfaces as well as the 
asymptotic geometry of the building blocks we describe below. We require the flexibility condition 
to guarantee the existence of a family of graphs that are smooth in two modifications we denote 
d, i. The parameter d will assign a vector to each vertex that is the sum of each unit direction 
of an edge or ray times its associated parameter f (2.6). The parameter £ will vary the length of 
each edge in the graph F. Notice that changing lengths may result in a change of direction for the 
edges and thus £ can influence d. We will be interested in graphs F that vary smoothly in d, £ for 
sufficiently small d, £. 

We will immerse an initial surface in R 3 based on the parameters that come from the graph F 
and from two additional parameters d, £• We call d the unbalancing parameter as it describes the 
deviation from balancing of a surface based on F and d (with £ = 0). An observation of Kusner 
regarding a homological invariant on CMC surfaces implies that CMC surfaces will be balanced, 
see [21, 22], and thus we presume F has d = 0. We call C the dislocation parameter as it will induce 
a dislocation of the surface from the structure induced by the graph. F and C together determine 
£ (4.16). 

Building blocks and the initial surface. The building blocks will be one of two types. The 

first type corresponds to an § 2 with geodesic disks removed and each block of this type will belong 

to a standard region on the initial surface. The second type of building block corresponds to an 

immersion of a piece of a cylinder that is a Delaunay surface except near the boundary. We designate 

the parameter of the family of Delaunay immersions by r and presume r G (— oo, 0) U (0, 1/4]. We 

let Y T : R x S 1 — > R 3 , see (3.5), denote the immersion of a Delaunay surface with parameter r. The 

asymptotic geometry of Delaunay surfaces is well understood; as r — > 0, the surface converges to a 

string of § 2 joined at antipodal points. Two types of geometric limits exist, where the type of limit 

depends on the sign of the Gauss curvature K. On a symmetric region where K > 0, as r — > the 

region converges to the round sphere. Rescaling by r -1 about the symmetric region where K < 

produces the standard catenoid - a complete minimal annulus in R 3 . We refer to both of these 

types of regions as standard or almost spherical regions. To justify this designation on a negatively 

\A\ 2 

curved region, consider the conformal metric h = —^~g, where g is the induced metric. Notice that 

2 



h is invariant under dilations of the surface and recall that for a minimal surface h is precisely 
z/*<7§2. In fact, in the metric h, the two geometric regions described are isometric. The convergence 
of the positively curved region to the round sphere thus justifies the name given to both regions. 

We designate a neighborhood of the region K = as a transition or neck region. In a second 
conformal metric, \ ( see (4-43)), each transition region is conformal to a flat [0,£] x S 1 where 
£ « — logr. In Figure 1 we identify the standard and transition regions on a Delaunay immersion. 
An extended standard region includes one standard region and the two adjacent transition regions. 

Extended Extended 




Figure 1. We label the regions of interest on a Delaunay piece. Notice both the 
positively curved and negatively curved standard regions look like S 2 in the h metric. 
Each shaded region becomes a long cylinder in the conformally flat x metric. 

We incorporate into the second building block a dislocation, designated by the parameter £. If 
one wishes to transit smoothly between a sphere and a Delaunay surface, the asymptotic geometry 
of a Delaunay immersion implies there exists a natural positioning for the two surfaces. The 
dislocation repositions the center of the sphere so the transit is no longer natural. In fact, to 
understand precisely the error terms, we apply two separate transitions. The first transition is 
applied to a Delaunay piece and a sphere placed in their natural positions, the second transition is 
from the original sphere to an off center sphere placed in a position designated by the dislocation 
parameter (Definition 3.20). 

The positioning of the pieces depends on the modified graph T(d, I) and the parameter £. At each 
vertex, we position a sphere with disks removed and at each edge or ray we position a Delaunay 
building block with axis depending on £. To guarantee the immersion is well defined, we must 
identify the neighborhood of each end of a Delaunay piece with the appropriate annular domain on 
an adjacent sphere. The initial surface has CMC identically one except on certain domains of each 
Delaunay building block, and the error is determined by the parameters r, £. 

The PDE and linearized equation. For a surface M C M 3 and a sufficiently small smooth 
function / : M — > R, it is well known that for Mf := {x + /(x)]Vm(x)|x G M} where Nm is a 
section of the normal bundle to M, one has 

H f = H M + (& + \A\ 2 )f + Q f . 

Here Hm,Hj represent the mean curvature of the surfaces M,Mf respectively, A and \A\ 2 are 
the Laplacian and norm squared of the second fundamental form in the induced metric, and Qf 
represents quadratic and higher terms in / and its derivatives. To obtain a CMC surface with mean 
curvature 1, we need to solve 

Cf := (A + \A\ 2 )f = l-H M -Q f . 

As in [6, 7, 9, 11-14, 16-19, 29], we solve the linearized problem on various regions of the surface in 
one of the conformal metrics h, x- The metric x se ts a natural scale by behaving - up to a small 
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reparameterization - like the flat metric on the cylinder for each Delaunay piece. In our particular 
setup, we determine 1 — Hm is small in the \ metric and thus we expect / should have small norm. 
The scale induced by x implies that the term Qf should not dominate. 

Our first goal is to prove Proposition 5.37, which provides a global solution to a modification of 
the original linearized problem Cf = 1 — Hm ■= E. Obstructions to solving the linearized problem 
correspond to the existence of kernel for the operator C. Additional technicalities arise because 
the process allows for small modifications to our surface which may induce small changes to the 
spectrum of the operator C. For this reason, we must identify the space of eigenfunctions for the 
linearized operator with small eigenvalue, a space we identify as the approximate kernel. On each 
standard region, we determine these potential obstructions to solving the linearized problem by 
taking advantage of the well understood geometric limit of each standard region in the h metric. 
Following the strategy of [6, 14, 20], we solve the linearized problem modulo the extended substitute 
kernel /C, a space of functions with properties outlined in Lemmas 5.23, 5.28. We introduce the 
extended substitute kernel to satisfy two criteria. On each extended standard region, we introduce 
a three dimensional space of functions in K, so that for each E there exists w £ /C such that E + w is 
1? orthogonal to the approximate kernel. This allows us to solve the modified semi-local problem 
on each extended standard region. We proceed to solve the global linearized problem by applying a 
global partition of unity to the inhomogeneous term. Since each transition region, in the x metric, 
is conformal to a long, flat cylinder, we solve the Dirichlet problem on each transition region with no 
modification to the inhomogeneous term. We then add these solutions to the semi-local solutions on 
their common regions of support to obtain a global one. To guarantee our global solution possesses 
good estimates, we must modify our semi-local solutions by prescribing the low harmonics on the 
boundary of an adjacent transition region to guarantee fast exponential decay on that transition 
region. This modification represents the second use of the extended substitute kernel. On each 
non-central standard region, the introduction of substitute kernel for the first use is also sufficient 
for the second. On each central standard region, we define a three dimensional space of functions 
on the core for the first use. For second use, we define a three dimensional space of functions at 
each attachment that modifies the solution along each of the attached standard regions. 



Prescribing extended substitute kernel and the geometric principle. At this stage, we 
have obtained a global solution to a modification of the linearized problem, thus solving one problem 
but introducing another. To correct this, we perturb the surface in such a way that we induce 
prescribed changes to the linearized equation. These changes take the form either of perturbations 
of the solution or perturbations of the initial surface. In both cases, our goal is to modify the mean 
curvature in such a way as to prescribe exactly the extended substitute kernel that was introduced 
in the previous step. Of course we cannot match the modification exactly, but we show the resulting 
error is small enough to close the argument. Small perturbations are accomplished by modifying 
the solution; larger perturbations require we modify the surface itself. At this step we rely on the 
flexibility of T which allows us to define a smooth family of immersions C M 3 which depend 
on T as well as d, C- 

We prescribe extended substitute kernel via the Geometric Principle [6, 14, 15, 18, 20]. On each 
non-central standard region, we use a balancing formula to prescribe substitute kernel. Balancing 
can be studied at the linearized level where it reduces to Green's second identity; thus the modifi- 
cation on these regions is done on the level of the function. On each central standard region, the 
parameters d, £ modify the surface in fundamentally different ways, although in both cases they 
control dislocations of the surface. As in [10] we introduce the parameter d to prescribe an element 
of the substitute kernel. Geometrically, it corresponds to an unbalancing of the graph V that results 
from repositioning the attached edges or rays at a vertex and changing each Delaunay parameter. 
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We introduce the parameter £ to prescribe the element of the extended substitute kernel that cor- 
responds to each attachment. As mentioned previously, at each attachment £ provides a vector 
that repositions the central sphere relative to each attached Delaunay piece. On the level of the 
function, the dislocation amounts to a linear combination of the Killing fields corresponding to the 
coordinates of the normal vector to the surface. We describe this relationship precisely in Lemma 
B.4. 

Prescribing the extended substitute kernel on the central spheres has some novel features because 
of the lack of symmetry in the construction. We prescribe the substitute kernel and extended 
substitute kernel in two separate steps and show that the error introduced in each of these steps is 
small enough to be absorbed. For more detail, see Proposition 6.9. 

Finally, using the geometric principle and careful estimates on the non-linear modification that 
we describe in Appendix A, we use Schauder's fixed point theorem to complete the proof of our 
main theorem (Theorem 7.1). 

Outline of the Paper. The paper consists of seven sections and two appendices. In Section 2, 
we describe finite graphs and their properties. We then provide examples of graphs that produce 
embedded surfaces. Section 3 includes a thorough description of the building blocks we use as 
well as estimates on the norm of the inhomogeneous term of the linearized equation. In Section 
4 we describe the immersion of an initial surface based on all of the parameters from T as 

well as d,C- I n Section 5 we develop the linear theory we need to solve the linearized problem 
Ch4> '■= (A/j + IA/jI 2 )^ = E. Finally, in Section 6 we describe how to modify the function and the 
surface to prescribe extended substitute kernel. We prove the main theorem in Section 7. Appendix 
A contains technical estimates corresponding to modifications of the initial surface. In Appendix B 
we define a function that describes the dislocation of the surface locally as a graph over the original 
surface. Additionally, we provide proofs for the estimates needed in the proof of the main theorem. 

Notation and Conventions. We set the convention that the mean curvature H := Kl + K2 where 
ki,K2 are the two principle curvatures for the surface. Thus, a unit sphere has mean curvature 
identically one. 

Throughout this paper we make extensive use of cut-off functions, and we adopt the following 
notation: Let \& : R — > [0, 1] be a smooth function such that 

(1) * is non-decreasing 

(2) ^ = 1 on [1, oo) and * = on (-oo, -1] 

(3) ^ — 1/2 is an odd function. 

For a,b £R with a ^ b, let ip[a, 6] : R — >■ [0, 1] be defined by ip[a, b] = ^ o L af) where L af) : R ->• R 
is a linear function with L(a) = —3,L(b) = 3. Then V'l )^] has the following properties: 

(1) *l>[a,b] is weakly monotone. 

(2) ip[a, b] = 1 on a neighborhood of b and ip[a, b] = on a neighborhood of a. 

(3) 4>[a,b]+ip[b,a] = 1 on R. 

Finally, it will be convenient for us to use weighted Holder norms and to that end define 

\\ip : C fc ' Q (tt, g,f)\\ : = sup -±-\\ip : C k ' a (B x D Q,g)\\ 

where Q, is a domain in a Riemannian manifold (M, g), f is a weight function, and B x is a geodesic 
ball centered at x with radius equal to the minimum of 1 and half the injectivity radius in the 
metric g. 
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2. Finite Graphs 



Our construction technique relies on properties of a finite graph that both determines the pa- 
rameters of various building blocks and describes the relationship between the building blocks. We 
begin by defining a more general class of finite graphs and then proceed to describe specific prop- 
erties we will need. The geometry of the building blocks will impose certain restrictions on lengths 
of edges in the graph, and the need for deformations will imply we are interested in graphs that 
have some flexibility. We impose additional conditions when we want to guarantee embeddedness 
of the surface in the main theorem. 

Definition 2.1. Let V denote a finite graph, by which we mean a collection 
{V(T),E(T),R(T),f} such that 

(1) V(r) is a finite collection of vertices, placed inM. . 

(2) E(T) is a finite collection of edges in M 3 , each with two endpoints in V(r). 

(3) R(T) is a finite collection of rays in M 3 ; each with one endpoint in V(r). 

(4) f : E(T) U R(T) -»• M\{0} is a function. 

We define two graphs as isomorphic if there exists a one-to-one correspondence between the 
vertices, edges, and rays such that corresponding rays and edges emanate from the corresponding 
vertices. 

Definition 2.2. Let E p denote the collection of edges and rays that have p G V(r) as an endpoint. 
We then have 

U peV{r) E p = E(T)UR(T). 

Also, for ease of notation we define the set 

(2.3) A(T) C V{T) x (E(F) U R(T)) 

such that [p,e] G -A(r) if e G E p . 

Note we choose the letter "A" here so that one thinks naturally of an attachment. 

Definition 2.4. For each e G E p , let v ejP denote a unit vector pointing away from p and parallel 
to e. For any p G V(r), let 

(2.5) V p = {v e , p |e G E p }. 
We let d T : V(T) -> M 3 such that 

(2.6) d r (p) = E f ( e Kp 

e&E v 

measures the deviation from balancing at the vertex p. 

Definition 2.7. If dr(p) = for all p G V(r), we say T is a balanced graph. 

To each edge of the graph we will associate a Delaunay piece with parameter rf (e) where r > 
is a global parameter that we choose close to zero so that rf(e) is small enough to meet all needs. 
(See Section 3 for a description of the Delaunay pieces.) Since the period of Delaunay surfaces 
tends to 2 as the parameter r tends to zero, we are interested in graphs with edge lengths close to 
multiples of 2. 

Definition 2.8. Let T be a finite, balanced graph. We say T is a central graph if each edge in 
E(F) has even integer length. Define I : E(T) — > Z + such that for e G E(T), 21(e) equals the length 
of e. 

Finally, we note the conditions needed to guarantee embeddedness. 
Definition 2.9. We say T is pre-embedded if it is a central graph with f : E(T)(J R(T) — > R + and 
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(1) For all p G V(T) and all ei / ej G E p , Z(v eiiP , v ejiP ) > tt/3, where Z(x, y) measures the 
angle between the two vectors x, y. 

(2) There exists e > such that for all e,e' G -E(T) U i?(r) that do not share any common 
endpoints, the Euclidean distance between e, e' is greater than 2 + e. 

(3) For any two rays e G E p , e' G E p > , 1 — v 6)P • v e /y > e. 

The conditions imposed on a pre-embedded T are the weakest we can currently allow to guarantee 
embeddedness of the CMC surface. The need for an e > in the second item comes from the fact 
that the maximum radius of an embedded Delaunay surface is on the order 1 — t + 0(t 2 ) but we 
allow for the edges to move with order Cr where C can be quite large. The first angle condition 
does not require an added e because the change in the period for r small (on the order — rlogr) 
dominates both the radius change and the changes we allow via unbalancing and dislocation (on 
the order r). The final condition is necessary because two parallel rays pointing into the same 
half-plane may cross in the process of unbalancing. 

Deforming the Graphs. In the construction process, we deform the central graph T via the 
parameters d, £. We need this deformation to be continuous and therefore are interested in graphs 
r that allow such continuous deformations. We describe how to deform a central graph based on 
two functions d, £. The function d will correspond to d/r and the function £ will depend on £, r 
and 1(e). Finally the order of i (—rlogr) will be much smaller than the order of d (order 1). 

Definition 2.10. Let L(T) be a space of functions from E(T) to M. with norm defined such that for 
I G L(T), 

\£(e)\ 

l = max 



eeE(r) /(e) 

7/r,r are isomorphic graphs with corresponding edges e,e', where T is a central graph, we define 
£r G L(T) such that 

2l{e)+2£ T {e) 

is the length of the edge e' . 

Definition 2.11. For a central graph V we let D(T) denote a space of functions from V(T) to 
M 3 . Let T,r be isomorphic graphs. Then D(T), D(T) can be identified in the natural way so that 
dip G D (r). Moreover we define, for any d G D (F), 

\\d\\ D = max \d(p)\. 
P ev(r) 

We now define the types of graphs we will use. 

Definition 2.12. Let V be a central graph. We say T is flexible if there exists e > such that for 
any £ : E(T) — > [— e,e] and d : V(T) — > {v G M 3 ||v| < e}, there exists a graph T(d,£) isomorphic to 
r and smoothly dependent on (d, £) such that: 

(1) Iff = T(d, 0) then df = d. 

(2) Iff = T(d,£) then £ f = £. 

Any flexible T will have an associated family of graphs .F(r). The existence of this family follows 
immediately from the definitions. 

Definition 2.13. We define a family of graphs, .F(r), to be a collection of graphs containing a 
central graph T, parameterized by a ball B(T) C D(F) x L(r), in a small neighborhood of the origin: 

(%)r = r(o,o). 

(2) F(d,£) is isomorphic to T(0,0) and depends smoothly on (d,£). 
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(3) Fori = r(d,0), d f = d. 

(4) ForY = T{d,£), t f = £. 

(5) Ife,e' are corresponding edges or rays for T(d,0),T(d,£) then 
( 2 - 14 ) ^o)( e ) = f W e ')- 

We remark that while (2.14) implies that the function f r ^ Q ^ equals the function f r ^^ on 

corresponding edges, the two graphs r(o?, 0) and T(d, £) may have different balancing at each vertex. 
That is, we do not expect that d r ,^ Jp) = d r ^ f ^(p') for corresponding p,p'. The change in 
unbalancing must be allowed as t will change the lengths of the edges and thus possibly the position 
of edges and rays about a vertex. The next lemma makes precise by how much the two unbalancing 
parameters can disagree. 

Lemma 2.15. Consider T(d,£) G J~(T) such that (d,£) are in a small neighborhood of the origin. 
Then there exists C = C(T) > depending only on the graph T such that 

(2-16) \\dw)-dmi)\\D<C(T)\\£\\ L . 
Proof. Notice that for a fixed p G V(r), 

where p^ G V(T(d, 0)),e^ G E p , and pi G V(T(d,£)), e? G E pt are corresponding vertices and edges 
or rays. Then the smooth dependence on £ implies that in a small neighborhood of the origin one 
immediately has 

l v e d -,p d ~ - v e/lW | < C\\£\\ L . 
As the graph is finite, the value f (e^) is absorbed into C(T) and the result follows. □ 

We now choose a frame associated to each edge in the graph T and use this frame to determine 
a frame on each edge for any graph in J 7 (T). 

Definition 2.17. For e 6 E(T) U R(T) we choose once and for all an ordered orthonormal frame 
Fr[e] = {vi[e], V2[e], V3[e]} such that Vi[e] is parallel to e. Moreover, for e G R(T), Vi[e] = v eiP . 

Because every edge is associated with two vertices, the choice of frame determines a direction 
along e G E(T) that allows is to distinguish between the two vertices. For convenience, we set the 
following notation. 

Definition 2.18. Fore G E(T) let p^\e\ be the endpoints of e such that 

v e,p+[e] = v lN = - v e,p-[e]- 

For e G R(F), let p + [e] be the endpoint of e. 
Also, let o e) p G {—1, 1} be defined by 

&e,p ■= v e , p • vi[e]. 

We define a general rotation that will be used repeatedly throughout the paper. 

Definition 2.19. Given two unit vectors x, y G M 3 such that Z(x, y) < it, let R[x, y] denote the 
unique rotation such that 

• I/xxy/0, let R[x, y] be the rotation that takes x to y by choosing the smallest rotation 
from x to y about the axis made by x x y. 

• 7/x = y, take R[x, y] to be the identity. 
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Finally, we use the rotation defined above to describe an orthonormal frame on any graph 
in the family T(T). Notice that as Z(v e , p ,v e 'y) < ir for any [p, e] £ ^4(r) and corresponding 
[p',e'] G A(T(d, £)), the rotation can be defined. 

Definition 2.20. For (d,£) € B(T) and e G E(T) U R(T) we define an orthonormal frame 
F r ^^[e] = {vi[e; d, £], V2[e; d, £], vs[e; d, £]} uniquely by requiring the following: 

(1) ~vi[e;d, £] = v e /y where e',p' are the edge (or ray) and vertex in the graph T(d,£) corre- 
sponding to e,p + [e] on V. 

(2) v, [e; d, £] = R[vi [e] , vi [e; d, £]] (v< [e]) fori = 2,3. 

Remark 2.21. F r 0^[e] depends smoothly on d,£. 

Examples of Graphs. As we wish to highlight the novel surfaces produced by our construction, 
we specifically mention graphs that produce embedded CMC surfaces. 

Our first examples are graphs that produce surfaces in Mo k where k G {3, . . . 12}; we do not 
need to presume any symmetry (though when k = 12 embeddedness is only guaranteed when the 
rays are placed symmetrically). We begin with T having one vertex and k rays emanating from 
that central vertex. The vertex position provides three free parameters and there are three free 
parameters for each of k — 1 rays. The balancing condition implies that this will entirely fix the A;-th 
ray both in magnitude and direction. Modulo rigid motions, there are 3k — 6 continuous parameters 
for each of these constructions. 
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Figure 2. Examples of pre-embedded T that produce surfaces in Mo,4, Mo$. 

We also consider graphs with a fixed number of edges and rays in Mo,k with restrictions on k 
coming from the number of edges in the graph. See Figure 2 for a few examples. In the first setup, 
we consider two vertices, one edge, and two rays emanating from each endpoint of the edge. The 
graph is not presumed to be planar, though we draw it as such for convenience. Choosing the 
position of the first vertex allows for three free parameters; with a fixed edge length, the position 
of the remaining vertex has two free parameters. We are free to choose f on the edge. Moreover, 
at each vertex we are free to choose one ray direction and f , resulting in six more free parameters. 
These choices and the balancing condition entirely fix the direction and weights of the remaining 
two rays. Modulo rigid motion, there are 6 = 3-4 — 6 free parameters and one discrete parameter. 
We provide another genus zero example in the second drawing in this same figure. We prescribe 
three edge lengths and position four vertices to accommodate these edges. The positioning of the 
vertices gives us 9 free parameters. Notice the edges need not be coplanar. To the interior vertices 
we add one ray and to the outer vertices we add two rays. If we prescribe f on each edge, the 
remaining 6 degrees of freedom come from the position and weight of two rays - one ray at each of 
the vertices that contain two rays. Modulo rigid motions we have 12 = 3 • 6 — 6 degrees of freedom. 
We note that each of these examples obviously demonstrate the required flexibility for performing 
a construction. We also note the pre-embedded conditions can be easily satisfied in each case. 

We can also construct an example of a finite, balanced, flexible, pre-embedded T with one edge 
and up to 22 rays - up to 11 rays at each vertex. With 11 rays at each vertex, embeddedness is 
guaranteed by positioning the rays and edge symmetrically (pointing to the center of the faces of 
a regular dodecahedron) and assigning the same parameter f to each edge and ray. We are free to 
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rotate one of the ray configurations about the edge and thus can guarantee no two rays are both 
parallel and pointing into the same half-plane. Moreover, for an edge of sufficient length, we can 
presume all the rays emanating from one vertex are of distance greater than 2 from rays emanating 
from the other vertex. 

By this method, we can easily produce infinitely many pre-embedded, flexible, finite central 
graphs that allow us to construct an embedded surface in A^o.fc- Generically, we may choose e 
edges and position e + 1 vertices to accommodate these edges. At each interior vertex, we can 
position between one and nine rays and at each final vertex we can position between two and ten 
rays. The limitation on the number of rays allows us to easily determine f on each edge and on the 
free rays so that the graph is balanced and pre-embedded. If we presume the necessary symmetry 
in the positioning of edges and rays, we can position 11 rays at the outer vertices and 10 rays 
at the interior vertices. We may have to increase the edge lengths in certain cases to guarantee 
pre-embeddedness. 

We can prescribe genus in the simplest way by introducing an equilateral triangle into the graph. 
For embedded surfaces, we require the equilateral condition. In Figure 3 we provide examples of 
central graphs that produce embedded CMC surfaces with genus 1, 2, and 3. For genus 1, we 
begin with an equilateral triangle which assigns to each edge the same length l(e\). Designate the 
parameter % > at e«. Then the position and weight of each ray is predetermined and we see the 
graph has three continuous and one discrete parameter. For genus 2, we begin with two equilateral 
triangles, sharing a common edge. We will position the minimum number of rays needed to satisfy 
the balancing condition presuming all weights are positive. We determine the edge length l(e\) and 
choose a weight fj > for each ej. Notice this provides us with five continuous parameters, but 
in fact we have one more. In this case, T need not be a planar graph. Thus, the angle between 
the planes containing the two triangles provides the final continuous parameter. With these six 
parameters, the balancing condition determines the position and weight of each of the four rays. 
We can construct an embedded genus 3 CMC surface in two ways. We use either a central graph 
T containing three equilateral triangles and at least five rays, or we can use a tetrahedral structure 
made out of equilateral triangles and at least four rays. We can also add more rays at each vertex 
which will increase the number of continuous parameters by three. In the case of three triangles and 
five rays, the graph T has one discrete parameter (edge length) and nine continuous parameters. 
Seven continuous parameters come from assigning a weight to each edge and the remaining two 
come by assigning an angle between the adjacent triangles. For the tetrahedron example, there 
exist six continuous parameters by assigning f j > for each edge ej. When we position one ray 
at each vertex, this completely determines the surface. When we position additional rays at each 
vertex, the number of continuous parameters goes up by three. Notice we can use the tetrahedral 




Figure 3. Graphs that produce surfaces in ^1,3,^2,4,-^3,4- 
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structure to build more non-planar finite graphs with v vertices and genus 2v — 5. From these 
graphs, we construct embedded surfaces in M2v-5,v 

Lemma 2.22. Given v > 4, there exists an embedded £ G M-2 V -5,v 

Proof. We proceed by induction. For v = 4, we choose a single tetrahedron made of equilateral 
triangles with edge length /. We put the same weight f > on each edge; the symmetry of the 
construction guarantees the rays point out radially in the direction connecting the barycenter (of 
the solid tetrahedron) with the endpoint of the ray. Therefore, no two rays are parallel and the 
graph is pre-embedded. Flexibility is easily see as any small d can be achieved by varying the ray 
weight and direction and changes to I can be accomplished by sliding vertices along the various 
edges or extensions of edges. 

At each stage of the process, we add one vertex and three edges, increasing the genus by two. 
Suppose for v vertices we have a flexible, pre-embedded, central graph T v . A quick calculation 
shows the genus of the embedded surface we construct is 2v — 5. We add one vertex v\ and three 
edges ei,e2,e3 to T v . We presume that each of the weights on the edges of T v remain fixed; the 
addition of the three edges provides three free parameters at this step. Notice that each edge ei 
has an endpoint at a vertex on r^, and an endpoint on v\. Thus, the choice of T{ ^> at each 
will influence the direction of four total rays. The pre-embedded conditions the rays must satisfy 
are the non-parallel condition and the distance condition. Thus, we have an open set of choices for 
the three fj and only finitely many conditions to satisfy. Therefore, we can easily determine fj > 
such that T^+i is pre-embedded. The flexibility also follows easily as again we can prescribe any 
small d by varying the ray direction and weight. Moreover, by adding a single tetrahedron to an 
already flexible graph, we are still free to smoothly vary t by sliding vertices along the particular 
edge of interest. 

Thus, for any v, we produce a pre-embedded, flexible, central graph T with v vertices and 3v — 6 
edges. The graph has 3v — 6 free parameters, coming from prescribing the f parameter to each 
edge. Moreover, the resulting embedded CMC surface has genus 2v — 5. □ 

3. The Building Blocks 

The surface we construct is built out of appropriately glued pieces of spheres and perturbed 
Delaunay surfaces. The positioning of these pieces and the parameter that describes the Delaunay 
immersion associated to each edge or ray rely on the graph T and the two parameters d, £. However, 
one can determine appropriate building blocks without any initial reference to the graph T. To 
highlight this fact, we first develop immersions of the building blocks that depend upon more general 
parameters but not on any graph. In Section 4 we explain how to use these general immersions, 
given a flexible, central graph V and parameters d, C to construct a complete surface in R 3 on which 
we prove the main theorem. 

Spherical Pieces. We first describe immersions of spheres with geodesic disks removed. In the 
initial immersion, these spherical pieces will be positioned at the vertices of the perturbed graph. 
Moreover, the centers of the geodesic disks will correspond to directions from which edges and 
rays of this graph emanate. As we are currently avoiding any reference to the graph, we consider 
immersions where geodesic disks are removed based solely on a collection of vectors. 
For a fixed, small S > 0, choose a such that 

(3.1) tanh(a + 1) = cos<5/8. 

Then for the sphere embedding we describe in Definition 3.12, Yo({ a + l} xS 1 ) C S 2 is the boundary 
of a geodesic disk of radius 5/8 centered at (1, 0, 0). 
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Consider a collection of unit vectors V = {xi, . . . , x ra } such that Z(xj,Xj) > 35 for i / j. For 
such collections of vectors, we will be particularly interested in the set 

(3.2) S 2 V = S 2 \ ]J Ef /8 {m) 

i=l,...,n 

where Df (x) denotes a geodesic disk on the unit sphere with radius r, centered at x. Notice 
the condition Z(xj,Xj) > 35 for i ^ j implies the distance between the removed geodesic disks is 
greater than 25. 

Proposition 3.3. Let V = {xi, . . . , x n } be a collection of unit vectors in M 3 such that Z(xj, Xj) > 
35 for i ^ j. If V = {yi, . . . y n } such that Z(xj, yj) < 5/16 for i = 1, . . . , n, then there exists a 
family of diffeomorphisms Y[V, V] : S v — > Sy, C R 3 , smoothly dependent on V , with the following 
property: 

Fori = l,...,n and x e Df /4 (xi)\Df /8 (xi) , 

R[x i ,y i ](x) = Y[V,V](x). 

For ease of notation, let Y[V, V] denote the identity map on Sy. 

Proof. For V, V as described, consider a smooth cutoff function ipy : Sy — > R defined so that 
4>v(x) = 1 for all x G § 2 \(UjD^ 3 (xj)), ipv(x) = 1 in a neighborhood of each <9D^) 3 (xj), ipv(x) = 
in a neighborhood of each D|^ 4 (xj) and smoothly transits between these values on each annulus 
centered at x,-. Define 



Y[V,V](x) :-- 



Y[V, V](x) if x e S 2 \ U Xi£ y £>f/ 2 (x*) 

^y(x)Y[V,V](x) + (l-^(x))R[xi, yi ](x) if x G Df /2 (xi)\Df /A (x t ), 

R[xi, yi ](x) iixeDf^XDf^) 



The main observation on smoothness comes by noting that the matrix R[xj,yj] is determined 
completely from smooth functions on Xj • yj and from the components of Xj x yj . Since these values 
are smooth in yj, the proposition holds. □ 

In order to deal with a small annoyance that will manifest when we identify pieces of the abstract 
surface M, we must introduce a second diffeomorphism. The process involves post-composing the 
previous map with a diffeomorphism that slightly rotates the image of the annuli under Y[V, V'\. In 
this way we guarantee that the diffeomorphism sends one specified orthonormal frame to a second 
specified orthonormal frame. 

Proposition 3.4. Let V\ x Vi = {(x];, x 2 ), . . . , (x^ , x 2 )} be a collection of pairwise associated 
unit vectors in R 3 such that xj ■ x| = for i = l,...,n and Z(x,J,xj) > 35 for i / j. Let 
V{ x V 2 = {(yj, yf ),..., (y*, y 2 )} be a second collection of pairwise associated unit vectors with 
yj • y 2 = for i = 1, . . . , n such that Z(xj,yj) < 5/16 for i = 1, . . . ,n and j = 1, 2. There exists 
a family of diffeomorphisms, Y_[Vi x V2, V[ x V£\ : S Vi —> Sy, C M 3 , smoothly dependent on V{, V^, 
with the following property: 
Fori = l,...,nandxe Df /4 (xj)\Df /8 (xj), 

Y[Vi x V 2 ,V{ x V&{x) = R\y[VuV{\(y*),y?] o Rfx^y, 1 ]^). 

Proof. The existence of such a diffeomorphism follows by using the same cutoff function and logic 
from the previous proposition. □ 

Observe that since Y[Vi, V/](x 2 ) is orthogonal to yj, the second rotation applied to each annulus 
fixes the vector yj for each i. 
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Delaunay Surfaces. While descriptions of Delaunay surfaces are readily available in the literature, 
our construction will require an understanding of our choice of conformal immersion into M 3 . To 
that end, we present, albeit in an abbreviated form, the conformal map and fundamental geometric 
quantities related to it. 

For r € (0, 1/4] define the immersion Y T : R x S 1 ->• M 3 

(3.5) Y T (t,9) = (k(t),r(t)cos0,r(t)am9) 

where 

r(t) = ^re w{ t\ 

fc'(t) = 2y/rr(t) cosh w(t) 
fc(0) = 0, 
with 

,„ fi x / (w') 2 + 4rcosh 2 u; = 1 

^' Dj \ w(0) > 0, w'(0) = 0. 

As a result, for any r € (0, 1/4], we have Gauss map, metric, and second fundamental form 

(3.7) u T (t,9) = (—w', 2y/r cosh w cos 9, 2>/r cosh w sin 9) , 

(3.8) g T = Te 2w (dt 2 + d9 2 ), 

(3.9) A T = -2Te w (smhw dt 2 + coshw d9 2 ). 

Notice for r < we replace each coshu; with a sinhw; and each r with |r|. 
With this definition, we see H = 1 and 

(3.10) \A T \ 2 = 2 + 2e~ 4w , K T = l-e~ 4w . 

Note the sign on the curvature of the surface corresponds precisely with the sign of w. 

By the nature of the equation and based on its initial conditions, w is periodic with period we 
designate \P T . Moreover, w has even symmetry about t = and odd symmetry about t = P T . It 
has a maximum at t = and a minimum at t = 2P T . The periodicity of w is preserved in the image 
surface, and we let 1 + p T = k(2P T ). Thus, the period of the image surface is 2 + 2p T . 



Lemma 3.11. A Delaunay surface with parameter r , as described above, is rotationally symmetric 
about the x\ axis and is embedded for r G (0, \]. Let r max ^ r mm denote the largest and smallest 
radii of the circles in the X2, X3 plane. Then 



Proof. The rotational symmetry and embeddedness follow immediately from the definition of the 
immersion and fact that k' > for r G (0, 1/4]. The ODE for w and initial conditions imply that 

for r > we have w(0) = arccosh (277) an< ^ w (?P r ) = — arccosh (277) ■ Using the log formulation 
for arccosh , 

r T (0) = v^e^W = (l + VT^) = \{^ + v 7 !^) = 1 + O(r), 

r T (2P T ) = jre w{2P ^ = (2y/r) (l + ^/l^y 1 = 2r (2 - 2r + 0(r 2 )) -1 = r + 0(r 2 ). 

Note when r < we replace arccosh by arcsinh and each r by |r|. Finally, any instance of — 4r 
becomes +4|r|. □ 

We will need a reference for an embedding of § 2 \{(±1, 0, 0)} from the cylinder and so state that 
here. 
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Definition 3.12. Let Y : R x S 1 ->• M 3 where 

Yo(t,9) = (tanhi, sechicosfl, sechtsinfl) . 
T/ien, c/o = sech 2 t(dt 2 + d6» 2 ) and |A)| 2 = 2. 

The construction relies on the fact that certain regions of Delaunay immersions possess well 
understood geometric limits as r — > 0. We solve the linearized problem with respect to a conformal 
metric h, which behaves on these regions much like the pull back of the Gauss map. In fact, u*g 
provides an isometry between the regions [—6, b] x S 1 and [2P T - b, 2P T + b] x S 1 via the map 
Y T (t,9) — > Y T (2P T — t,6). (Here b is a large constant, fixed in (5.3).) Therefore, it suffices to 
understand the asymptotics of the immersion of [—6, b] x S 1 . 

In [10], the regions and their geometric limits are described in some detail. The next lemma will 
be stated without proof. The interested reader should consult Lemmas 2.1 and 2.2 in Appendix A 
of [10] for the details. 

Lemma 3.13. Let ] — > R be the function whose graph, rotated about the x\ axis, gives 

Y T ([-2P T , 2P T ] x § 2 ) . Then as r 

p T —5-0, X T — > 1. 

Let ro(xi) : [—1, 1] — > R be defined by ro(xi) = \J\ — x\. Given e > 0, there exists r e > such that 
if < \t\ < T t , then r T restricted to [—1 + e, 1 — e] depends smoothly on r and 

\\r - r T : C k ([-1 + e,l - e})\\ < C(e,k)\r\. 

Moreover, we have the following period limits as r — > 0: 

(3.14) i im _^.*r = i ; i im Pr 



r^o — logr dr t-^o — rlogr 

We have the following corollary comparing the metric on the sphere and the Delaunay immersion. 

Corollary 3.15. For e G (0, 1), k G Z + , there exists T e > such that for all < |r| < r e , 

(3.16) \\ gT - go: c k ([-l + e,l-e]x8 l ),g T )\\<C(e,k)\r\. 

For a fixed large b (choose for example the largest b such that tanhft = 1 — e), Lemma 3.13 
expresses the limit, as r — > 0, of the immersions of the regions [AP T n — b,4P T n + b] xS 1 for each 
n£Z. Regions of the form [2(2n — 1)P T — b, 2{2n — 1)P T + b] x S 1 are isometric to these regions in the 
metric v*g T under the mapping Y T (t, 9) — > Y T (2P T — t,9). On regions in between, one cannot appeal 
to natural geometric limiting behavior. Instead, we understand the behavior of these portions of 
the cylinder in the flat metric dt 2 + d9 2 . We determine the limiting length of such a cylindrical 
piece in this metric. 

Lemma 3.17. There exists b » 1 such that 

IP — lb 

(3.18) lim ^ M = 1. 

r^o — log r 

Proof. For e' > 0, (3.16) provides b sufficiently large (for any < |r| < tv small) such that 
-w' T 6 (1 — e', 1] on [b, P T ] x S 1 . Thus, we can determine, as w(P T ) = 0, 

fPr 

(P T - b)(l - e') < - w' T dt = w T {b) < P T - b. 
Jb 

For b fixed, recall that \r (b) - r T (b)\ = \A - tanh 2 b - ^r~e w ^ < C{b)r. Thus, 

log (Vl - tanh 2 b - C(6)r) 2w (b) log (\/l - tanh 2 6 + C(6)r) 

1 + = < P ' < 1 + = 

— log r — log r — log r 
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Letting r — > proves the result. 



□ 



Delaunay Building Blocks. We now describe a general immersion of an appropriately perturbed 
Delaunay piece. A verbal description of the immersion follows Definition 3.20. Throughout this 
subsection we presume a few technical conditions on parameters of interest, which will always be 
stated in the hypotheses. These conditions are satisfied for the setup in Section 4 and are necessary 
to attain the desired estimates. Throughout this subsection, let < r" C 1 be small enough to 
ensure the immersions described below are smooth, well-defined immersions and let C denote a 
large universal constant independent of r" . 

We first define four cutoff functions that will be useful. 

Definition 3.19. Let ii dislocation^ A gluing^ ■ (t,6) G [a,4P T l - a] x S 1 ->• R 3 be cutoff functions 
such that: 

• *Pdislocation+ = ^[a + 2, a + 1], 

• ^dislocation- = ^[4P T / - (a + 2), AP T l - (a + 1)], 

• i>gluing+ = ^[« + 3, a + 4] , 

• ^ gluing- = ~ (a + 3), AP T l - (a + 4)]. 

With these cutoff functions we define the immersion. Observe the smallness conditions on r, 
defined at the outset are to ensure the immersion is a smooth, well-defined immersion. 

Definition 3.20. Given r G M\{0},a G R+,l G Z+,C ± G M 3 with < (C^ < Ct",0 < |r| < 
t", 1 <C a, we define two smooth immersions Y e d ge [T, I, C + , C~] : [a,4P r Z — a] x S 1 — > M 3 and 
Y ray [T, C + ] '■ [«> oo) x S 1 R 3 such that, for X = (t, 0), 

Y ed ge[T, I, C + ,C]{%) =^dislocation+(t){ Y o(x) + C + ) + (1 ~ *l> dislocation (*)) (1 ~ W (*))^0 (») 

+ 1 Pgluing+ if) ' ^ 'gluing- {£) ' Y T (x) 

+ (1 - ^dislocation- (*))(! ~ f'sl-r (*)) y ~ ( x ) + ^dislocation- (*) (^o" (») + C~) 

^ray[T,C + ](^) =, 0di S / OC ati O n+ (*) (^0 (a?) + C + ) + (1 ~ ^dislocation+ (*) ) ( 1 ~ *Pgluing+ (t))Y (x) 
+ i>gluing+(t) " ^r(^) 

w/iere F ~(x) = F (* - 4P r Z, 0) + (2 + 2p T ) Zei. 

To aid the reader, we describe the geometry of the Y e d ge immersion in some detail. For t G 
[a, a + 1], the image is a geodesic annulus sitting on a unit sphere centered at £ + . The annulus is 
centered at £ + + ei with inner radius 5/8. When t G [a + 1 , a + 2] , the immersion smoothly transits 
between an annular region on the dislocated sphere and an annular region on a unit sphere centered 
at the origin. For t G [a + 2, a + 3], the immersion remains on the unit sphere centered at the origin, 
while for t G [a + 3, a + 4], the immersion smoothly transits between this sphere and a Delaunay 
piece with parameter r. The same procedure happens toward the other end. First, the Delaunay 
piece transits back to a unit sphere centered at (2 + 2p T ) lei. This position represents the location 
of the end of a Delaunay piece with parameter r and I periods, with initial end on the {x± = 0} 
plane and axis on ei. Finally, this sphere transits to a unit sphere centered at £~ + (2 + 2p T ) le±, 
a dislocation of £~ from the previously described sphere. 

Of course, the Y ray immersion has the same behavior as Y e d ge near the origin. The only difference 
is that the Delaunay immersion continues out to infinity and there is no transiting back to a sphere. 

In determining the initial immersion, the unbalancing parameter d will induce changes in the 
r parameter associated to each edge or ray. As we stress here the independence of the building 
blocks on the background graph, we consider a diffeomorphism of a piece of the cylinder that will 
account for length changes based on a change of parameter associated with unbalancing. Again, 
this diffeomorphism can be determined from parameters that are independent of any graph V. 
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Definition 3.21. Given t,t' such that |r| < t" and 

(3.22) J G (1 - Cr", 1 + Cr") 

/or some (7 > 0, we define two families of diffeomorphisms 

Y ed ge 

[r, t'] : [a, AP T l - a] x S 1 ->■ [a, 4P T /Z - a] x S 1 

and 

Y ray [T, r'] : [a, oo) x S 1 ->• [a, oo) x S 1 
which depend smoothly on t,t', such that 

Y edge [r, r'](t, 0)=i;[a + 5,a + 4](t) • (t, 0) 

+ ^[a + 4, a + 5](t) • Y>[4P T Z - (a + 4), 4P T / - (a + 5)](t) • 



+ ^[4P r Z - (a + 5), 4P T / - (a + 4)](t) • (i + 4Z(P T , - P T ), 9) . 

and 

Y ra y[r, r'](i, 6) = i/j[a + 5, a + 4](t) • (t, 0) + V[a + 4, a + 5](t) ( 



Again, the conditions on r, r' in the definition above will always hold for immersions of interest 
to us. Moreover, the following lemma makes clear that by presuming these conditions we can state 
C fc -norm bounds in Proposition 3.28 without involving 

Lemma 3.23. 



i p " 



Pr 



-logr 



(3.24) 

Proof. The estimate follows immediately from (3.22) and (3.18). □ 

Definition 3.25. We define two new immersions, smooth in all parameters, Y_ edge [r, r', I, £ + , £~] : 
[a, 4P T l - a] x S 1 -»• M 3 and Y ray [r, r', C + ] ■ [a, oo) x S 1 -»• R 3 suc/i toai 

(3.26) Y edge [ T , t', I, C+, C"](i, 0) = Y edge [r', I, C + , C] ° ^e[r, r'](t, 0) 

and 

(3-27) y ray [r, r', C + ](t, d) = Y ray [r', C + ] o tayfr r'](t, 0). 

Proposition 3.28. Lei g := (y e d se )*<7R3 or g := (Y_ ray )*g^:i as the situation dictates. For a fixed, 
large constant b > a + 5, 

\\Y edge [T,r',lX + X-} ~ Y : C k ((a,b) x S\<7)|| < C(M)(|C + I + \r\) 

\\Y edge [r,r',lX + ,C}- Y - : C k ((4P T l - b, 4P T l - a) xS\g)\\ < C(k,b)(\C\ + \r\) 

and 

\\Y ray [r,r'X + ] ~ Y : C fe ((a,6) x S\g)\\ < C(k,b)(\C + \ + |r|). 

Proof. The statement follows immediately from the definition of the immersions, (3.16), and the 
conditions on r, r' which allow us to invoke (3.24) as needed. □ 

Proposition 3.29. The image surface Y_ edge [r, r', /, £ + , C~]([°> 4P T Z — a] x S 1 ) has mean curvature 
identically one except when t G [a + 1, a + 2] U [a + 3, a + 4] U [AP T l - (a + 2),4P T l - (a + 1)] U [4P r i - 
(a + 4),4P T Z- (a + 3)]. 

The image surface Y_ ray [T, r', £ + ]([a, oo) x S 1 ) /tas mean curvature identically one except when 
t G [a + l,a + 2] U [a + 3, a + 4]. 

16 



Proof. The statement is obviously true for the immersions Y edge ,Y ray based on their definition. 
Moreover, Y edge ,Y ray do not change the mean curvature of the immersions Y edge , Y ra y d 

Definition 3.30. Let H[edge] and H[ray] represent the mean curvature of the surfaces immersed 
by Y_ edge and Y_ ray respectively. Let H error [edge] = H[edge] — 1 and H error [ray] = H[ray] — 1. 
Further, let 

H dislocation\'\ + Hgluing['] = E error [-] 

where 

supp (H dislocation [edge}) C ([a + 1, a + 2] U [AP T l - (a + 2), 4P T Z - (a + 1)]) x S 1 , 
SUpp (-H'disJocationt'""!/]) C ([a + 1, a + 2]) x S 1 , 

supp (H gluing [edge\) C ([a + 3, a + 4] U [4P T i - (a + 4), 4P T l - (a + 3)]) x S 1 , 

supp (-H"^ U m 9 [ray]) C [a + 3, a + 4] xS 1 , 
where supp (/) represents the support of the function f. 

With this definition, Proposition 3.28 provides the following corollary. 
Corollary 3.32. Let g := (Y_ e dgeY '9r 3 or 9 := (Y-ray)*9R 3 as the situation dictates. Then: 

(1) \\H glumg [edge]-.C k ([aAPrl-a]x^),g)\\<C(a,k)\r\. 

(2) \\H gluing [ray] : C k ([a, oo) x S 1 ),^ < C{a,k)\r\. 

(3) \\H dlslocatlon [edge] : C k ([a, AP T l - a] x S 1 ), g)\ \ < C(a,k) (|C + | + |C"I)- 

(4) \\H dls iocaUon[ray\ : C k ([a, oo) x S\<7)|| < C(a, k)\t + \. 

4. Construction of the Initial Surface 

Given a flexible, central graph T and parameters d, £, we determine an immersion into R 3 by 
appropriately positioning the building blocks described in Section 3. Because our construction relies 
on the geometric limit of certain pieces of a Delaunay surface, we need the parameter associated 
to each e G E(T) U R(T) to be small. To that end, we choose a parameter < r <C 1 and define 
e G E(T)UR(T), 

r (e) := rf (e). 

The argument requires that all |ro(e)| be sufficiently small so that we may appeal to various geomet- 
ric estimates. Thus r will depend on maxf, but not on d, £, or the structure of T. The finiteness of 
r implies there exists a fixed constant C such that t/C < |ro(e)| < Cr. We follow the convention 
that parameters with a "hat" are on scale one and those with no hat are on scale r. 

The Abstract Surface M. Given a flexible, central graph V with rescaled functions ro(e), we 
define an abstract surface M that will serve as the domain for the immersion into M. 3 . M depends 
only on T and is independent of the parameters d, C- The construction of M proceeds as follows. 
To each vertex p G ^(r), we associate M[p], a two-sphere with geodesic disks removed. The centers 
of these geodesic disks correspond exactly with the directions of the vectors v eiP G V p . The radius 
of the disks depends upon a S we choose based on the graph. To each e G E(T) UR(T) we associate 
a piece of 1 x S 1 with left boundary at {a} x S 1 where a depends on 5 as in (3.1). If e is an edge, 
the length of the piece is determined based on the functions ro(e) and 1(e). Finally, we define the 
abstract surface M by appropriately identifying a geodesic annulus on M\p] centered at v ejP with 
a neighborhood of the appropriate boundary of the cylindrical piece M[e\. 

For the graph T, fix a small < 5 <^ 1 such that for each p G V(F) and e, e' G E p , Z(e, e') > 3d. 
We now fix once and for all the constant a so that 



(4.1) 



tanh(a + 1) = cos(5/8). 
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Definition 4.2. For p G V(r) define 

(4.3) M[p] = S 2 Vp 

where this notation follows that of (3.2) and V p is defined in 2.4- As the enumeration along a 
cylinder corresponding to an edge will vary depending on the direction, we simplify notation by 
defining 

(4.4) P[e] := 4P To(e) /(e). 
For e G E(T), let 

M[e] = [o,P[e]-a] xS 1 

while for e G R(F), let 

M[e] = [a,oo) x S 1 . 

To make the proper identification between M[e] and M[p], for e G E p , we define a rotation that 
takes the standard frame in M 3 to the frame Fr[e] associated to e G £ p . 

Definition 4.5. For e G E(T) U R(T), let R[e] : M 3 -»• R 3 denote the rotation such that 

R[e](ei) =Vi[e] 

/or i = 1,2,3, where here the Vj[e] re/er to i/ie ordered orthonormal frame chosen in Definition 
2.17. 

With this rotation in hand, we define the abstract surface M. 
Definition 4.6. Lei 

m'= ( n ) n( n m n 

\pev(r) / \eeE(T)uR(r) 
and let 

(4.7) M = M' I ~ 

where we make the following identifications: 

For [p, e] G A(T) with p = p + [e] and x G M[e] n ([a, a + 1] xS 1 ), 

x~ (R[e]oy (i))nM[p]. 

For [p,e] G A(T) with p = p~[e] and x = (t,6) G M[e] D ([P[e] - (a + l),P[e] - a] x S 1 ), 

x~ (R[e]oy (t-P[e],S))nM[p]. 

Standard and Transition Regions. Following the strategy in [6, 9, 10, 12-18, 20] we carefully 
define various regions on M. We require slightly more burdensome notation than in previously 
cited work as our setup lacks any symmetry. For each p G V(r), the region S[p] is a central 
standard region [6, 18] or a central almost spherical region in the terminology of [9, 10, 12-17]. The 
geometric limits of certain regions of the immersion Y T : R x S 1 — > M. 3 motivates our identification of 
various regions on each M[e]. Recall that 21(e) denotes the length of each edge on the initial graph 
r. Thus, each M[e] will have 21(e) — 1 geometrically well understood regions and 21(e) regions 
connecting them. This leads us to define the following sets. 

Definition 4.8. We let [p, e, n] G S(F) where 

S(T) :={[p, e, n]\e G E(T), [p + [e],e] G A(T),n G {1, 2, . . . , /(e)}} 

U {[p, e, n]\e G E(T), \p~ [e], e] G A(T),n G {1, 2, . . . , /(e) - 1}} 
U {[p, e, n]\e G R(T), [p, e] G A(T),n G N}. 
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Further, let [p, e, n'] G iV(r) w/iere 

iV(r) :={b, e, n']|e G £7(r), [^[e], e] G A(r), n' € {1, 2, ... , /(e)}} 
U {b,e,n']|e G P(r), [p,e] G A(T),ri G N}. 
Finally, we let \p,e,n"] G A^ + (r) where 

N + (T) :={\p,e,n"]\e G £(r), [p^e] G A(T),n" G {0, 1, . . . , 1(e) - 1}} 
U {[p, e, n"]\e G R(T), \p, e] G A(T), n" G N U {0}}. 

With this notation, each S\p, e, n] C M will correspond to a standard region or almost spherical 
region. Each A[p, e, n'] will correspond to a transition or nec/; region. Notice that <S(r) C N(T). 
In fact, iV(r)\5(r) = {[p _ [e],e,Z(e)] : e G -E(r)}. We choose this notation so that the set 5(r) 
enumerates every standard region exactly once. Moreover, the enumeration of the standard regions 
is such that it increases along M[e] as one moves further away from the nearest boundary. For 
e G E(T), the middle standard region on M[e] bears the label ( S'[p + [e], e, 1(e)]. Each S\p, e,n] is an 
extended standard region and contains both the standard region and the two adjacent transition 
regions. The S\p] are central extended standard regions and contain all adjacent transition regions, 
where adjacency is determined by e G E p . The C + \p,e,n"],C~[p,e,n'] represent the meridian 
circles on the boundary between a standard and transition region. C[p, e, n'] is the center meridian 
circle of A[p, e, n'\. 

Recall that a is determined by 5 in the definition of M and controls the size of the removed 
geodesic disks. The constant b determines the size of each standard and transition region. We 
use x, y in subscripts to modify the size of the regions and the boundary circles. For example, 
S\p] C S x [p] while S x [p] C %]■ 

Definition 4.9. Let p G V(r). For S(T), N(T), N + (T) as outlined above, we define the following 
regions on M. 

(1) S x \p] := M\p] U {eb=p+[e]} (M [e] n [a, b + x] x S 1 ) 

U{e\p=p-[e]} (M[e] n [P[e] - (b + x),P[e]-a] x S 1 ) 

(2) S x [p] := M[p] U {eb=p+[e]} (M[e] D [a, 2P To(e) - (b + x)] x S 1 ) 

U{e|p=p-[e]} (M[e\ n [P[e] - (2P ro(e) - (6 + x)),P[e] - a] x S 1 ) 
^ 5 x [p + [e], e, n] := M[e] D [2nP To(e) - (6 + x), 2nP To(e) + (6 + x)] x S 1 
^ 5 x [p"[e],e,n] := M[e] n [P[e] - (2nP To(e) + (6 + x)),P[e] - (2nP ro(e) - (6 + x))] x S 1 
^ 5 x [p+[e], e, n] := M[e] D [(2n - 2)P To(e) + (6 + x), (2n + 2)P To(e) - (b + x)] x S 1 
(6) S x \p-[e],e, n] := M[e] n [P[e] - ((2n + 2)P T()(e) - (6 + x)), P[e] - ((2n-2)P To(e) + (b + x))} x S 1 
^ A^ip+fe], e, n'\ := M[e] n [(2n' - 2)P To(e) + (6 + x), 2n'P ro(e) - (6 + y)] x S 1 
^ A x , y \p-[e],e, n'\ := M[e] n [P[e] - (2n'P To(e) - (6 + y)),P[e] - ((2n ; - 2)P T()(e) + (b + x))] x S 1 
^ C+[p+[e], e, n"] := M[e] n {2n"P To(e) + (6 + x)} x S 1 
fiflj C+[p-[e], e, n"] := M[e] D {P[e] - (2n"P T()(e) + (5 + x))} x S 1 
^i; C-[p+[e],e, n'] := M[e] D {(2n'P ro(e) - (b + x)} x S 1 
(12) C-]p-[e],e, n>] := M[e] n {P[e] - ((2n'P To(e) - (6 + x))} x S 1 
C[p + [e],e,n'] := M[e] D {(2n' - l)P To(e) } x S 1 
C[p-[e],e,n'] := M[e] D {P[e] - ((2n' - l)P T()(e) )} x S 1 
i/ere 6 > a + 5 is chosen in (5.3) and is independent of r i»Me < x,y < P T — b where positivity 
of P T — b is guaranteed by the smallness of t. We set the convention to drop the subscript x when 
x = 0; i.e. S\p] = S x \p\. Moreover, we denote K XyX = A x . For the sake of clarity, we point out: 

( V S x \p\ = S x [p] |J ee Ep A x\p,e, 1] 
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Non-central 
\ Standard 
Region 



Extended Standard Region 
of a Central Sphere 



C-\p-[e],e,l(e)} 



\s\p+[e],e,l(e)} 



\C-[p+[e},e,l(e)} 



A[p-[e],e,Z(e)] 




A[p+[e],e,/(e) 



Enumeration at the middle of an edge e 

Figure 4. Two schematic renderings of M. The top one is near a vertex p with 
\E p \ =3 and the bottom one at the standard region associated to center of an edge 
e. Note that standard regions appear spherical and transition regions appear cylin- 
drical. 



(2) S x \p,e,n] = A x \p,e,n] U S x \p,e,n] U A x \p, e, n + 1] for all n £ {1, . . .,l(e) - 1} 

(3) S x \p, e, 1(e)] = A x [p+[e], e, /(e)] U S x \p, e, /(e)] U A x \p-[e],e, /(e)] 

(4) C+\p,e,0} = S x \p]nA x \p,e,l] 

(5) C+\p, e, n"\ = S x \p, e, n"\ n A x [p, e, n" + 1] /or n" € {1, . . . , /(e) - 1}_ 

(6) C~ [p, e, n] = £3; [p, e, n] n A x [p, e, n] except when n = 1(e) and p = p~ [e] in which case 

C^1e],e,I(e)] = S I ^[e],e,!(e)]nA I [p-[e],e,I(e)] 

(7) dS x \p]=U ejeEp C+\p,e j ,0] 

(8) dS x \p, e, n] = C x \p, e, n] U C x [p, e, n] for n G {1, . . . , /(e) — 1} 
f9j aS x [p,e,Z( e )] = C-[p + N,e,/(e)]UC-[p-[e],e,Z(e)] 

dA x , y {p, e, n'] = C+ [p, e, n' - 1] U C" [p, e, n'] 

The graph r(cZ, £). In this subsection, we explain how to modify the initial graph T by prescribed 
parameters d, £. The modification by d is immediate while the modification by £ induces £. As 
T is flexible, there exists a family T(T) such that for (d,£) in an e neighborhood of the origin, 
F(d, t) G T(T) varies smoothly in d, I. 
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Because we scaled the function f by r, we do the same for d and thus we are interested in 
d G D(F) such that 

(4.10) \\d\\ D <er. 

Here e is taken from Definition 2.12. For any such F(d, 0), let : E(F) U R(F) — > R denote the 
rescaled f r(j>0) . 

We now determine the function £ : E(T) — > R that will rely - for each e - on /(e), Td{e), and two 
vectors Cb + [ e L e L C[p~ [ e L e ] G M 3 . The parameter £ describes a dislocation of each Delaunay piece 
from its central sphere and thus we describe the full parameter as a map from each attachment. 

Definition 4.11. Let C ■ A(T) ->• R 3 such that 



(4.12) 



C\P, e] 



i=i 



For i/ie graph T, we define C £ I^ 4 ? /inife dimensional vector space with values in R, indexed 
over i = 1, 2, 3 and \p, e] G ^4(r). Zei 

IICII = . max Kb, e] I . 

\p,e]eA(F) 

As we will see, the norm of £ can be quite large compared to the norm of d. Throughout the 
paper, we allow 

(4.13) IICH <Qt 

where C is a large, universal constant that is independent of r. 

To each edge in the graph we associate a Delaunay building block with length determined by 
the parameter r^(e) and the number of periods /(e). Let / : E(F) — > R such that 

(4-14) T(e):={2 + 2p Td{e) )l(e). 

Thus, a Delaunay piece with /(e) periods and parameter r^(e) will have axial length equal to /(e). 
Definition 4.15. For e G E(F) we define the function £ : E(F) ->• R by 
(4.16) 2£(e) := \c[ P + [e], e] - (cb~M, e] + (7(e), 0, 0)) | - 2/(e). 



Y + C + 




Y ~ + C" 



Figure 5. In the figure, we let ( + ,( correspond to Cb + [ e ]i e ]iCb [ e L e ] respec- 
tively. Also, notice that Y ~ is defined so that its center is at (2 + 2p Td ^)l(e). 

The definition of £{e) amounts to the following steps (see Figure 5). First, we position a segment 
of length /(e) so that it sits on the positive xi-axis with one end fixed at the origin. Then we 
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dislocate the two ends of this segment corresponding to C[p + [ e ]> e ] an d C[p - [ e ]> e ] where Cb + [ e ]> e ] 
is the dislocation from the origin. We then measure the length of the segment connecting these two 
points. Finally, we compare that length with the length of the edge e in the graph T. 

Lemma 4.17. For £ G L(T) with £(e) as defined in (4.16) and r sufficiently small, 
\W)\\ L := SI <CCr-Cr d (e) log r d (e) < -Ct log r. 



Proof. Notice 



-c||cil+p Td(e) <^|<c||cil+p Td(e ). 



Thus, applying (3.14) and the uniform bound on ||£|| immediately implies the result. □ 
We now highlight a few flexibility and asymptotic results we will use throughout the proof 

Lemma 4.18. For r sufficiently small and \ \d\\D < £T, ||C|| < Qj we have the following: 

(1) Ife,e' are corresponding edges or rays for T,T(d,0) then there exists some constant C such 
that 

(4.19) € (1 — Ct, 1 + Ct) , and Z(e, e') < Get 



(4.20) 



(2) and 

Ct 



P r d (e') 



< 



logr 



P r (e) 

(3) For T(d,0),T(d,£) and corresponding edges e,e', 

(4.21) Z(e,e') < -Crlogr, and |K(d,o) ~ dv{d/)\\D < -Ct 2 logr. 

Proof. The proof follows from the flexibility of T, an appropriate modification of Lemma 2.15, and 
Lemma 4.17. The estimate (4.20) follows immediately from (3.24) and the uniform bounds on the 
ratio of |Td(e)| and r. □ 

Remark 4.22. Notice that the finiteness of the graph T, the definition of To(e), and the condition 
(4.19) implies that for any \\d\\D < £T, there exists C such that \/C < J^^li < (j xhis gives us 
the freedom to replace any bounds in (^(e)! 1 * 11 by (Ct) ±1 , reducing notation and bookkeeping. 

The Smooth Initial Surface. We describe an immersion of the initial surface M which depends 
upon the initial graph V and the parameters d, £. The immersion amounts to an appropriate 
positioning of the building blocks described in Section 3. As is evident from the notation, the 
building blocks themselves depend upon the parameters and the graph. 

The positioning of the building blocks will depend in part upon a rotation associated to each 
edge that takes an orthonormal frame of the edge connecting 1(e) + C[f - [ e ]> e ] an d C[p + N> e ] to the 
orthonormal frame of the edge e' 6 E(T(d, £)) corresponding to e G E(F). 

Proposition 4.23. ForC, from Definition 4-11 and each e £ E(T) there exists a unique orthonormal 
frame F^[e] = {ei[e], e2[e], e3[e]}, depending smoothly on C, such that 

(1) ei[e] is the unit vector parallel to (c[p~ [e], e] + (1(e), 0, 0)^ — C[p + [ e ]> e ] such that e± [e] ■ e\ > 
0. 

(2) Fori = 2,3, e;[e] = R[ei,ei[e]](e;). 

(3) For any v G M 3 , |v - R[ei,ei[e]](v)| < C||C|| • |v|. 
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Proof. The first two items can be done by definition. To see the last item, observe that 

|ei x ei[e]| = |sin Z(e l5 ei[e])| := |sin0[e]| < C||C||. 

As |v - R[ei,ei[e]](v)| < \/2\v\yfl - cos0[e], the result follows. □ 

Definition 4.24. For e G R(T) we simply let ej[e] = ej. 

With this frame, we can describe the rigid motion that will position each building block associated 
to an edge or ray in R 3 . 

Definition 4.25. For each e G E(T) U R(T) with e' denote the corresponding edge or ray on the 
graph T(d,£), let R'[e] denote the rotation in M 3 such that R'[e](ej[e]) = Vi[e;d,£] fori = 1,2,3. Let 
T[e] denote the translation in R 3 such that T[e](R'[e](C[p+[e], e])) = p + [e']. Let U[e] = T[e] o R'[e]. 

Then for all a G R, 

(4.26) U[e] (C[p + [e},e\ + c iei [e\) = p + [e'} + Ci v t [e; d, I]. 

The building blocks positioned at vertices of the graph T(d, £) correspond to spherical pieces with 
geodesic disks removed. The rigid motion required for positioning these in M 3 is simply a translation, 
but the building blocks themselves are determined by a diffeomorphism that depends upon the 
initial frame Fr[e] and the new frame F^[e] for each e 6 E p . 

For each p G V(T), let {ei, . . . , e\ Ep \} be an ordering of the edges and rays in the set E p . Let 

(4-27) Vp 1 x V p 2 = {(v ei , p , v 2 [ei]), . . . , (v e|Bp| , p , v 2 [e| Ep |])}. 

This set of ordered pairs represents the direction the edge e emanates from p and the second element 
in the frame .Fr[ e ]- Notice that if p = p + [e] then the ordered pair is actually the first two elements 
of the frame Fr[e] and if p = p~[e] then the first element in each pair corresponds to — Vi[e] while 
the second element is unchanged. Recalling Definition 2.18, let 

V^' 1 = {<Te llP R'[ei](ei), • • • , ^^R'te^iKei)}; V£' 2 = {R'[e 1 ](e 2 ), . . . , R'[ e | £p |](e 2 )}. 

These sets represent the first two elements of the frame that will position the building block corre- 
sponding to an edge or ray emanating from p. One should note that in general R'[e]ej / Vj[e;cZ, £}. 
Figure 6 may help to highlight this fact. Thus, the geodesic disks removed under the diffeomor- 
phism do not correspond to the frame F r r d j\ [e] but to another frame that is dislocated from this 
one depending on the vectors Cb + [ e L e L C[p~ [ e ], e ]- With these tools in hand, we define the initial 
immersion. 

Definition 4.28. Let Y^,C : M — > M 3 be the smooth immersion defined by 
(4.29) 

' P' + YXVp 1 x T/ 2 , Vy x V p C \x) for x G M\p] 

Y dti {x) = \ U[ e ]oy e(i9e [ro(e),T d (e),i(e),C^[e],e],C[p1e], e ]](x) for x G M[e], e G E(F) 
{ U[e] o Y ray [T (e),T d (e),C\p, e]](x) for x G M[e],e G R(T) 

where here p' G V(T(d,£)) is the vertex corresponding to p G V(r). 
We denote 

(4.30) M dX = Y dX (M) 

and let H d £ : M — > R denote the mean curvature in the induced metric. 

Remark 4.31. When T is pre-embedded, the embeddedness of the image surface M d £ follows imme- 
diately from the smallness of r, the embeddedness of each of the building blocks and the smallness 
of d, C 
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Figure 6. A rough idea of the immersion of one edge. Notice that the transfor- 
mation U[e] sends the dislocated spheres to the vertices of the graph. The bold 
segment in the bottom picture corresponds to the positioning of the edge on the 
graph F(d,£). The Delauney piece has axis parallel to R'[e]ei, not parallel to the 
corresponding edge on the graph F(d,£). 



On each M[p], embeds a sphere with geodesic disks removed where the centers of these disks 
match with the directions (7 ep R'[e]ei for e 6 E p . The twisting diffeomorphism that describes the 
building blocks M[p] must be included so the immersion is well defined on M - in particular on 
the regions M [p] n M[e\. 

The immersion of each M[e] amounts to positioning a building block of a Delaunay piece from 
Section 3, where the parameters of the building block and its positioning depend on £ and F(d,£). 
The rigid motion of this piece positions the Delaunay portion of the immersion with axis parallel 
to R'[e]ei. Notice the Delaunay piece may not have axis parallel to the associated edge of the 
graph T(d,£). In fact, the piece will only have axis parallel to the associated edge of F(d,£) if 
C[p + [e], e] = C[p - [ e ]i e L which implies ei[e] = ei. 

Definition 4.32. For e G E(F) define H dislocation^], Hgiuing[ e ] '■ M — > R in the following manner: 

H MM ._ / H dA - 1 ifx G M[e] n (([a, a + 2] U [P[e] - (a + 2), P[e] - a]) x S 1 ) 

n dislocation^) .- | Q otherwise, 

n MM ._ j H d)i - 1 if x £ M[e] n (([a + 3, a + 5] U [P[e] - (a + 5), P[e] - (a + 3)]) x S 1 ) 
n 9luing[e\{X) .- | Q otherwise. 
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For e £ R(T) define 

_(H dX -l i fxeM[e}n([a,a + 2]xS 1 ) 
n dislocati<mie\{x) — | Q otherwise, 

tt H dX -l ifx€M[e]n([a + 3,a + 5] xS 1 ) 

nglvingWW — y Q otherwise. 

As an immediate consequence of the immersion and the definitions we have the following corollary. 

Corollary 4.33. All of the functions described above are smooth. Moreover we can define Hdisiocation '■- 
Y.eeE(T)vjR{T) H dislocation[e\ and H g i uing likewise. Then the smooth function H error := H d ^ - 1 : 
M — > R can be decomposed as 

Herror H dislocation ~t~ H g [ u ^ ng . 

Geometry. We introduce maps Y\p] and Y\p, e, n] on S\p], S\p, e, n], respectively, which one can 
interpret as the geometric limit of the normal of Y d ^ as r — > 0. A comparison between these 
maps and Y d ^ will accomplish two things. First, in the next subsection, we use these limits to 
compare two conformal metrics h, \ to the pull back by the Gauss map of the induced metric g. 
Second, using techniques from Appendix B of [10], these maps help us understand the kernel of 
the linearized operator on Y d £ by considering the linearized operator in the induced metric of Y 
(which corresponds to the metric on the round sphere) . 

Definition 4.34. Let Y[p] : S\p] -»• S 2 such that for x = (t, 6), 

{ Y[V p ,V p ](x) tfx€M[p], 

R[e] o Y (x) ifp = p + [e],xe M[e] n ([a,2P ro(e) - b] x S 1 ) , 

R[e]oy (t-P[e],fl) ifp = P~[e], 

x G M[e] n ([P[e] - (2P To(e) - b), P[e] - a] x S 1 ) . 

Let y [p, e, n] : S'fp, e, n] — >■ S 2 suc/i £/ta£ 



(4.35) Y\p]{x) 



The definition of V \p, e, n] when n is odd has the effect of reflecting the map when n is even over 
the {x\ = 0}-plane. As we wish to compare the normals of Y\p, e, n] and a rotated Y d £, this change 
is necessary as the normal for Y d £ when n is odd is a reflection of the normal when n is even. 

Proposition 4.37. Let p 6 V(r) and let p' be the corresponding vertex in the graph F(d,£). Then 

\\(Y d , c -p')-Y{p} :C\SM,{Y\p]T9&)\\<C{k, X yl 2 . 

For e e E(T) U R(T) 

||U[e]- 1 oJV diC (x)-y|p, e ,n](s):C fc (^b,e,n],(y|p,e,n]^ 

Proof. Recall by definition 

Y edge [T O (e),T d (e),l,C\p + [e],e],C\p-[e],e]](t 1 0) = 

Y e dge[Td{e),l, C [p + N , e] , C \p~ [e], e]] o y edse [T (e), r d (e)](t, (9). 

By (4.20), the reparameterization function will only impact derivatives on this function by multi- 
plicative factors of |= which are bounded by 1 — Therefore, we consider Y edge ,Y ray in place 
of Yedgei¥-ray throughout this proof as the bounds on the estimates we determine can absorb the 
multiplicative factors. 
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We consider first the behavior of Y\p] on the region M[e] n S x \p] for e G £? p . For e G -E(r), the 
C fe bounds are controlled by bounding 

(R'[e] oY edge [T d (e),l,C\p + [e\,e\X[p'[e],e}} - R'[e]oY ) + (R'[e]oY - R[ei, v e , p ] o Yb) • 

Recalling Proposition 3.28, the C k bounds of the first parenthesis are controlled by C(k, x)(\ \£\ \ + 
\Td(e)\). For the second parenthesis, observe 

|(R'[e] - R[ei, v e , p ])v| < C(\\d\\D + \\£\\l + ||CII)|v|. 
As (4.10), (4.13), and Lemma 4.17 imply 

l|d|b + ll% + ||CII < (Q. + e)r - Ct log r < r 1 / 2 , 
the result follows. (Similar arguments hold for Y ray .) The diffeomorphism on M\p] that describes 
the immersion depends smoothly on ||£|| and thus the C k bounds on M[p] follow immediately. 

For the second norm, first notice that for a general Delaunay immersion Y T , the evenness of w' 
about P T and (3.7) imply v T {t,6) = v T (2P T — t,9). Thus, it is enough to prove the inequality for 
any n even as the description of the function Y[p, e, n] for n odd accounts for the reflection. Recall 
that on S x \p,e,n], for e G E(F), 

U[ e ]- 1 oy di< (x) = y e(i9e [ro(e),r d (e),/,Cb + [e],e],C[p-N,e]](t,^. 
Using (3.16) to compare the immersions Y Td ^ and Yq on S x {p, e, n] for n even implies the result. □ 

Conformal Metrics. Following [6,10,14,18-20], we solve the linearized problem and prove the 
main theorem by appealing to conformally changed metrics. The metric h will be useful for under- 
standing the approximate kernel as every standard region will behave like an S 2 in this metric. We 
prove the main theorem in the % metric, a conformal dilation of the induced metric, which provides 
a natural scale for the immersion. 

We first define a cutoff function that will be helpful for a global description of the metric %. 

Definition 4.38. Let : M — > [0, 1] be a globally defined smooth function on M such that ip = 1 
for x G S[p]nikf[j)], ip := tp[a + 5,a + 4:](t) for (t,9) £ S\p]nM[e] andip is identically zero elsewhere. 

Thus, ip(x) is identically 1 on each central sphere S\p] except in a neighborhood of its boundary, 
where tjj smoothly transits to zero. Notice ip, 1 — if> form a subordinate partition of unity on M. 

Definition 4.39. We define the function p : U e6 £(nuii(r)-^[ e ] - ► K such that for each (t, 9) G M[e], 

p(t,e) = r Td(e) (t_,9)- 1 . 

Here we define 

(4.40) t = -p^t. 

Let p : M — > M be the globally defined smooth function such that 

(4.41) p(x) = Hx) lA ^ + (1 - 4(x))p(x). 

Definition 4.42. Let g = (^d,^)*ffR3. On M we define two global metrics h,x such that 

I A\2 

(4.43) h =-^T9-, X = P 2 9- 

We now list a number of metric relations that will be useful throughout the paper. 

Lemma 4.44. Assuming r is small enough, the following hold: 
(1) \\h:C k (S x [p],g)\\<C(k,x) 
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(2) \\h-(Y\p])*g S 2) : C k (S x \p], (Y\p])*g § 2)\\ < C{k,x)r^ 2 

(3) \\ X :C k (S x \p],h)\\<C(k,x) 

(4) \\h-(N d , c )*g:C k (S x \p,e,n},h)\\<C(k,x)T 

(5) \\h - (Y[p,e,n]rg S 2 : C k (S x [P,e,n],(Y\p,e,n])*g s2 )\\ < C(k,x)r 

(6) \\ X :C k (S x [p,e,n],h)\\<C(k,x) 

(7) \\p ±l :C k (S x \p],x)\\<C(k,x) 

(8) Wp^-.^iM^P^W^Cik) 

Proof. The definition of h and the uniform geometry of the region S x \p\ for r small implies (1); (2) 
follows immediately from Proposition 4.37 and the triangle inequality. Lemma 3.13, Proposition 
4.37, and the definition of the metrics h and % immediately imply (3). 
Using (3.10), on each S x \p, e,n], 

(4.45) fc=(W^)x. 

Further, h = 2\r d {e)\ cosh(2w(t))(dt 2 + d6 2 ) and x = df + dO 2 . Thus, (N dX )*g, h, and x all provide 
a natural isometry between positively curved standard regions and negatively curved standard 
regions on each M[e]. Therefore, it is enough to prove the next three items on a standard region 
with limit § 2 , that is, a standard region with Gauss curvature K > 0. Item (4) is proven in 
[10], Appendix A. We adapt the statement there to fit with our notation. For (5) we appeal to 
Proposition 4.37 and (4). Item (6) follows from Lemma 3.13 and (4.45). 

On S^fp], (7) follows from the definition of p, Proposition 4.37 and the uniform equivalence 
of g,h,x on each £x[p]- For (8) we define a new function and consider its C k bounds. For any 
x = (s,6) G M[e]\S'[p] let p(x) := j^fj^ for some fixed i G M[e]\5[p]. By definition, 

p( x ) = e -w(s)+w(t)_ 

We are interested, of course, in C k bounds on p{x) ±l for x G (i - 3/2, t + 3/2) x S 1 := D. Recall 
the bounds on p Td(e) imply that D contains a unit ball about the point (i, 9) £ M[e] with respect 
to the x metric. The definition of w implies \w'\ < 1 and thus 

\w(s) — w(t)\ = 

Thus, Hp* 1 : C°(D,x)\ \ < C. Now 

o 
-p{x) = -w'(s)p(x); gpP(x) = -w"(s)p{x) + (w'(s)) 2 p(x). 

The definition of w implies w"(s) = — 4\Td(e)\ cosh w(s) sinhw(s) = |Trf(e)|(e 2 "'^-^ — e~ 2w ^). Re- 
calling that C\Td(e)\ < r Td ( e ) < 1 + C\Td(e)\ on each M[e], ^(e) = \ T d{e)\ l l 2 e w , and w is odd 
about P Td ( e ) we observe that \w"\ < C on each M[e\. The uniform bound on \w'\ along with the 
previous analysis implies | |p ±1 : C 2 (D,x)\\ < C for some constant C. For any k G N, we can write 
^p(x) as a function of ^mw{x) and p{x) for m = 0, . . . , k — 2. This follows simply from repeated 
differentiation of the ODE of w and the fact that p is an exponential in the variable t. Thus, we 
obtain C k upper bounds for p on each M[e] via the C k ~ 2 upper bounds. Moreover, for derivatives 
of p -1 , we apply the same logic and note that derivatives will now include factors of /5 _fe_1 , which 
is bounded as a function of k. Putting all of this together with item (7) implies item (8). □ 

Finally, as we choose to solve the global linear problem in the metric x> we state precisely the 
estimates we have for the error on the mean curvature with respect to the metric x- 
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w'(t)dt 



< C\t- 



Proposition 4.46. \\p- 2 H gluing : C k (M, X )\\ < C(k)r. 
\\p- 2 H dlslocatlon : C k (M, X )\\ < C(k)\\C\\ < C(k)Cr. 

Proof. The bounds follow immediately from Proposition 3.32 and item (7) of Lemma 4.44 using 
the uniform equivalence of the metrics g, \ on the support of the functions involved. □ 

5. The linearized equation 

The goal of this section is to clearly state and prove Proposition 5.37. Recalling (4.43), we start 
by defining the operators 

C g := A g + \A\ 2 ; C h = j^C g ; C x = p' 2 C g . 

We choose to study the inhomogenous linearized equation in one of the following forms: 

C x 4> = E- C h 4> = T^ E - 

The linear equation on the transition regions. We first consider the linearized equation 
on A XtV \p, e, n'] - the transition regions. Throughout this section, we refer to each of these regions 
simply by A and their boundary components, C+ [p, e, n' — l] and C~ [p, e, n'] respectively as C + , C~ . 

On the transition regions, % takes the form dt 2 + dO 2 . That is, x is the flat metric on the 
reparameterized cylinder. Let x, x, x denote the t-coordinate distance from C + ,C~ and dA = 
C+UC". Then 

(5.1) (2n' — 2)P T0 ( e ) + b + x + x = t, 2n'P ro ( e ) — b — y — x = t, x_:= min{x, x}. 
Let i\ define the t-coordinate length of the cylinder parameterizing A. Then 

(5.2) £ A = 2P To(e) -2b-x-y. 

We consider weighted Sobolev spaces in the x metric where the weight function is with respect 
to x. This poses no problem as the weight function is a decaying exponential, and (4.20) implies 
that exponential decay with respect to t and t differ only by multiplicative factors of the form 
e T = l + 0{ T ). 

For each ei > 0, there exists b » 1 so that for A defined in terms of that b, 

(5.3) ||p- 2 |A| 2 :C' fc (A,x,e^)||<6 1 . 

This follows from the work done in Lemma 3.17, once we observe p~ 2 \A\ 2 = re 2w {2 + 2e~ 4w ) = 

At cosh(2w). Recall t = p d< ~ e) t. Each A with parameter To, equipped with the metric x, is isometric 

r o( e ) 

to [b, 2P TQ — b] x S 1 equipped with the metric dt 2 + d6 2 . Moreover Lemma 3.17 and the fact that 
t/C < to < t implies 

(5.4) |^A+logT|<C 

where C depends on e\ and on T. 

We now choose ei > appropriately small to satisfy the desired conditions going forward and b 
large enough to satisfy (5.3). For this fixed 6> 1, 

Proposition 5.5. The lowest eigenvalue for the Dirichlet problem for C x on A is > Ct~^ . 

Proof. The proof follows exactly as in [14], in particular refer to the proof of Lemma 2.26. Using 
the Rayleigh quotient to determine the lowest eigenvalue gives the lower bound described. □ 

Immediately, we get the following corollary. 

Corollary 5.6. (1) The Dirichlet problem for C x on A for given C 2,/3 Dirichlet data has a 
unique solution. 
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(2) For E G C ' 13 there is a unique ip G C 2 ' 13 such that C x ip = E on A and ip = on <9A. 
Moreover : C 2 ^(A, X )|| < C{f3)i\\\E : C°^(A,x)||. 

Definition 5.7. Let T-L^C] denote the finite dimensional space of spherical harmonics on the merid- 
ian sphere C that includes all of the functions in the first k-eigenspaces of A§i. 

Observe %o[C] is the space of constant functions and Hi[C] is spanned by {1, cos 6>, sin#}. In the 
next proposition and corollary, we examine the Dirichlet problem where we are allowed to modify 
the low harmonics on one boundary circle. This induces fast exponential decay for the norm of the 
solution to the Dirichlet problem on the transition region. 

Proposition 5.8. Given (3 G (0, 1) and 7 G (1, 2) there exists e\ > such that the following holds. 
There is a linear map K A : C°>P(A) C 2 '^(A) such that the following holds for E G C°^(A) and 
V := K A E: 

(1) C X V = E on A. 

(2) V\c+ G Hi [C + \ and vanishes on C~ . 

(3) \\V : C 2 'P(A, X ,e~r-)\\ < C \\E : C°^(A, X , e~^)\\. 

(4) TIa depends continuously on all of the parameters of the construction. 

The proposition still holds if the roles of C + and C~ are reversed and x is replaced by x. 

Proof. The first two items follow by standard linear theory. For the third and fourth, let TZa denote 
the linear map applied to the operator A x and suppose V\ is such that A x Vi = E, i.e. IZaE = V\. 
Then 

WC^-E: C°^(A, x ,e-^)|| < C(/3, 7 )ei||£ : C°^(A, x , e~^)\\. 

For ei sufficiently small, C x is a small perturbation of A x and thus an iteration implies the result. 

□ 

Corollary 5.9. Assuming ei is small enough in terms of 62 and (3 G (0,1), 7 G (1,2), there is a 
linear map: 

n d :{fe C 2 ' /3 (C+) : / is L 2 (C + ,d6 2 )- orthogonal to fti[C+]} -»• C 2,/3 (A) 

such that the following hold for f in the domain oflZg and V := TZgf ': 

(1) C X V = on A. 

(2) V\ c + — f G Hi [C + \ and V vanishes on C~ . 

(3) \\V\ C+ - f : C 2 >P(C+,d8 2 )\\ < e 2 \\f : C 2 ^(C+, d6 2 )\ \ . 

(4) \\V : C 2 ^(A, X , e -^)|| < C{M)\\f ■■ C 2 ^(C + ,d9 2 )\\. 

(5) IZg depends continuously on all of the parameters of the construction. 

The proposition still holds if the roles of C + and C~ are exchanged and x is replaced by x. 

Proof. Let 

n d :{fe C 2 ' /3 (C+) : / is L 2 (C+, d# 2 )-orthogonal to Wi[C + ]} -»• C 2 ^(A) 

define the linear map such that for / in the domain and V = IZgf one has 

(1) A X V = on A. 

(2) V\c+ = f and V vanishes on C~ . 

(3) \\V : C 2 ^(A, X ,e^)|| < C(f3,j)\\f : C 2 >?(C+, d6 2 )\ |. 

Such a map exists by standard linear theory. The corollary follows immediately by applying the 
previous proposition to the map 

K d f := Kgf - K A C x K d f. 

□ 
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We consider the behavior of solutions to the Dirichlet problem with prescribed boundary data 
in %i[C + ] U %i[C~} for the operator C x by comparing these solutions to known solutions for the 
operator A x . As we frequently reference and use these solutions, we provide a thorough definition 
and description. 

Definition 5.10. Fori = 1,2,3, let Vi[A, 0,1,0,2], Vi[A, ai, 0,2] denote the solutions to the Dirichlet 
problems on A given by 



with boundary data 



£ x Fi[A,ai 


a 2 ] = 0, 


A x ^[A,oi, 


02] = 


--0, 


Vi[A, ai,a 2 ] 


= Vi[A,ai 


a 2 ] = ai 


on 


C+, 


Vi[A,ai,a 2 ] 


= ^i[A,ai 


a 2 ] = a 2 


on 


c~, 


V 2 [A, ai,a 2 ] 


= ^2 [A, ai 


a 2 ] = ai cos 6 


on 


C+, 


V 2 [A,a 1 ,a 2 ] 


= V 2 [A, ai 


0-2] = a 2 cos# 


on 


c~, 


V 3 [A,a 1 ,a 2 \ 


= K 3 [A,ai 


a 2 ] = ai sin 


on 


C+, 


y 3 [A,ai,a 2 ] 


= V 3 [A,ai 


a 2 ] = a 2 sin 9 


on 


c~. 



As Vi[A, ai,a 2 ] and Vi[A, a\,a 2 \ are linear on ai,a 2 and the roles of C + ,C can be exchanged, it 
is enough to understand Vi[A, 1,0], Vi[A, 1,0]. We observe 
(5.11) ' 

VI A, 1,0 = — , V 2 A, 1,0 = . , , p , /p > cos6>, V3 A, 1,0 = . ^ ; °^ ; sing. 

sinh(P rcj(e) £ A /P ro(e) ) smh(P Td(e) ^ A /P To(e) ) 

Proposition 5.12. V"i[A, ai,a 2 ] is constant on each meridian and V 2 [A, a±, a 2 ], V^A, ai, a 2 ] are 
constant multiples of cos sin# respectively on each meridian. Moreover, for e\ sufficiently small 
in terms of €3 > 0, /3 G (0, 1),7 G (1,2), £/iere are constants Af, A± , A 2 , A3 such that 

(1) \Af - 1| < e 3 ,j4r < e 3^- 

f2j ||Vi[A,l,0]-VI[A,A+ Ai] iC^CA^.e-^ + r^-^JH < e 3 . 

(S; |A 2 -1| <e 3 ,JA 3 -l| <e 3 . 

^ \\V 2 [A,1,0] - V2[A,A 2 ,0] : C 2 ^(A, x , e~(^ + e^)|| < e 3 . 

f5j ||V 3 [A,1,0] - K3[A,A 3 ,0] : C 2 ^(A, x , e"^ 1 ^ + e^)|| < e 3 . 

Proof. The rotational invariance of the solutions follows from the rotational invariance of the op- 
erator C x . 

By (5.3), and the definition of Vi the decay takes the form of: 

||£ x Fi[A,l,0] : C^(A, X ,e-^ + £- l e~^)\\ < e 1 C(f3, 7 ) 

and 

I \C X V X [A, 0, 1] : C°>"(A, X , + r V^)! | < £l C(/3, 7 ). 

We now apply Proposition 5.8, but with the weaker estimates available here. Upon application, we 
determine the existence of V\[A, 1, 0], V1[A, 0, 1] G C 2,/3 (A) such that 

C X V X [A, 1, 0] = -C X V X [A, 1, 0]; C X V X [A, 0, 1] = -C X V X [A, 0, 1], 

both functions are constant on C + ,C~ by rotational invariance, and finally 

||Vi[A, 1,0] :C fl ^(A, X ,e-^ + r 1 e-^)||< e 1 C(P, 1 ) 

and 

|| V1[A, 0, 1] : C 2 ^(A, X , e-T* + rV^JH < eiC(/3, 7 ). 
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We choose At , A x so that 

V^A, 1, 0] = Fi[A, A+, A^} + A+Vi[A, 1, 0] + A^V^A, 0, 1] 

is achieved on dA = C + U C~ . Notice both sides of the equation are in the kernel of C x on A. 
The uniqueness of the Dirichlet solution on A implies equality holds on all of A. The estimates on 
At, Aj come simply from applying the given estimates on the appropriate meridians. 
For Vi[A, 1, 0] with i = 2, 3, we determine 

||^[A,1,0] :C^{A, X ,e-^ 1 >- + e- i )\\<e 1 C{P, 1 ). 
The rest of the proof follows the strategy outlined previously. 

□ 

Corollary 5.13. If f £ C 2,/3 (A) satisfies C x f = on A and f\ c - = 0, then 

(1) ||/ : C^(A, X , (* + l)/£)\\ < C(b,f3)\\f ■■ C^(C+,d8*)\\. 

(2) \\d t f : C°(AnC,x)\\ < C(b,0)\\f ■■ C^(C+,d0 2 )\\/\\ogT\. 

Proof. On C + , decompose / = f\ + ft + + f 1 - where /i = ai , /2 = «2 cos 0, f^ = 03 sin and f 1 - 
is in the domain of IZq defined in Corollary 5.9. Let V = T^df ± given by Corollary 5.9. Then using 
the estimates from Corollary 5.9 and Proposition 5.12 implies (1). Using this same decomposition 
and analysis, (2) follows by recalling (5.4). □ 

The approximate kernel on the standard regions. We now consider the nature of the ap- 
proximate kernel for the operator Ch on the standard regions. By approximate kernel, we mean the 
span of eigenfunctions of that have eigenvalues close to 0. Following the methods of Appendix 
B in [10], we understand the approximate kernel by comparing each standard region with the h 
metric to a round sphere with the standard metric. Because we do not impose symmetry, the span 
of the approximate kernel is three dimensional on each standard region. We determine a basis for 
the approximate kernel on each standard region using a function defined on that domain of M. 
The first functions defined below are normalizations of the components of the unit normal to the 
immersion Y^. The second are also based on the unit normal to the immersion, but are defined 
relative to a Delaunay building block with axis positioned on the xi-axis. 

Definition 5.14. Let ft: M R for i = 1, 2, 3 such that 

(r. 1 r.\ i ei ' Nd >C 1 at 

(5 - 15) f ^ : =\^N^m = ^- N ^ 

where Nd£ is the normal to the immersion Y^. 

Let fi : U e6 E(r)uH(r)-^[ e ] ~~ * K for i = 1, 2, 3 such that for x G M[e] 

(5-16) h{x) := llei . Nsi: v m ( = -«i ■ (R N N^x)) . 



Recall that the functions 



N&2 • e ?: 



|| ei -iV §2 :L 2 (S 2 )|| 

form an orthonormal basis for the kernel of the operator A§2 + 2. This observation and an un- 
derstanding of the geometric limiting objects (in the h metric) of the standard regions as r — > 
motivates the definition of the functions. The incorporation of a rotation in the definition of fi 
makes the choice of constants in Definition 5.27, (6.1) and the estimates they satisfy more imme- 
diately obvious. 



Proposition 5.17. For fixed h large, let e > 0. There exists r t such that for < r < r e , 
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(1) Ch acting on S\p,e,n], with vanishing Dirichlet conditions, has exactly three eigenvalues in 
[— e, e] and no other eigenvalues in [—1/2, 1/2]. There exist {/i[p, e, n], f2\p, e, n], /s[p, e, n]} 
i/tai denote an orthonormal basis of the approximate kernel where fi\p, e, n] G C°°(S'[p, e, n]). 
Moreover fi[p,e,n] depends continuously on the parameters of the construction and each 
fi\p, e, n] satisfies 

(5.18) \\fi\p,e,n] - h : C 2 ^(5 5 b, e, n], h)|| < e. 

fi?) £ft acting on S[p] has exactly three eigenvalues in [— e, e] and no other eigenvalues in 
[—1/2,1/2]. There exist {/i[p],/2[p],/3[p]} that denote an orthonormal basis of the ap- 
proximate kernel where fi\p] £ C°°(S\p\). Moreover, fi\p] depends continuously on the 
parameters of the construction and satisfies 

\\fi\p]-f i :C 2 >f ) (S 5 \p],h)\\<e. 

Before we proceed with the proof, we make the following relabeling of transition regions. 

Definition 5.19. For any standard region S\p,e,n], define A c / ose ,Aj ar such that 

Kiose ■= A[p,e,n] 

A _/A[p,e,n + l] if n< 1(e), 
A /--\ A[ p -[e],e,l(e)} if n = 1(e). 

Notice that the wording of this definition makes sense everywhere except on the middle of each 
edge. That is, except when n = /(e), the labeling of "close" and "far" correspond to relative 
distance to the nearest boundary circle on each M[e]. At the center of M[e], we keep the labeling 
but note that both transition regions are closer to a boundary than the standard region between 
them. 

We now prove the proposition. 

Proof. We prove both items in the proposition by taking advantage of the results of [10], Appen- 
dix B, which determine relationships between eigenfunctions and eigenvalues of two Riemannian 
manifolds that are shown to be close in some reasonable sense. Before proceeding, note that the 
inequality B.1.6 in [10] should read 

||^/||oo<2||/||oo. 

To simplify the language of the proof, we let S denote a general S\p, e, n]. Moreover, let 
C 'dose, C f ar denote the center meridian circle of A c i ose , Af ar respectively. Finally, let Y ,Y c i ose ,Yf ar 
denote the immersions defined in (4.36) for the particular [p, e, n] of interest and those that also 
share A c i ose , Af ar , respectively, on their extended standard regions. Thus, for 2 < n < 1(e), 
Ydose '■= Y\p,e,n — 1], and for 1 < n < 1(e) — 1, Yf ar := Y[p,e,n + 1]. Moreover, when n = 1, 
Ydose ■= Y\p] and when n = 1(e), Y far := Y\p~[e],e,l(e)]. 

We begin by considering S\ (^C close U Cf ar ^j . Let S + denote the connected component containing 
S and let A^ ose ,Aj ar denote the two remaining connected components where naturally we have 
^dose ^ A c ; ose . Let D c i ose , D f ar . denote the smallest geodesic disks in S 2 that contain, respectively, 
Yciose(A~ lose ) , Yf ar (Aj ar ). Consider two abstract copies of S 2 , {1} x S 2 ,{2} x S 2 and their disjoint 
union {1,2} x S 2 . 

We first determine a comparison Riemannian manifold (S, g) on which the linear operator of 
interest is well understood. In order to apply the results of [10], we must show there exist linear 
maps such that all of the assumptions of B.1.4 are satisfied. To that end, let 

S := ({1} x S 2 ) |J ({2} x (D close U D far )) 
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and define g on {1, 2} x § 2 by making g the restriction to g§2 on each copy of § 2 . Consider the 
operator A§2 on S 2 with the standard metric. It it well known that the lowest eigenvalues for 
this operator are 0, 2, 6, where the eigenvalue 2 has multiplicity 3. Therefore, the only eigenvalue 
of A§2 + 2 in the interval [—1,1] is zero, with multiplicity 3. Choosing b sufficiently large, thus 
ensuring the radii of D c i ose , Df ar are sufficiently small, we can guarantee the lowest eigenvalue of 
Ag on (D c i ose j ar , g) with Dirichlet data must be larger than 3. 

Define £ : S -> S such that Y(x) := (l,Y(x)) if x G S + and Y(x) := {2,Y doseJar (x)) if 
x G ^ciose far- We now construct two functions Fi,F2 that will satisfy the assumptions of B.1.4 
in [10], thus allowing us to appeal to Theorem B.2.3 in that same paper. First, let ip be a cutoff 
function on S such that ip = 1 on 5 and ip := ip[2d,d] o x on A c ; ose , A f ar . This provides a 
necessary logarithmic cutoff function that yields appropriate gradient bounds; see B.1.8 in [10]. 

For / G C§°(5), define o Y = ^ • / and for / G Cg°(£) let F 2 (/) = ^ ■ / o Y~\ For d 

sufficiently large, the assumptions of B.1.4 are easily checked and one can apply Theorem B.2.3 
with r = 2, n = 4, d = 1. That is, for the operator A§2, there is a gap between the first and second 
eigenvalues and between the fourth and fifth eigenvalues (counting multiplicity). 
Now recall that an orthonormal basis of the kernel of A§2 + 2 are the functions 

h > S2 \\e t .N s2 :L 2 (S 2 )||' 

i = 1,2,3, where e« are the standard unit basis vectors for M 3 and N S 2 is the outward normal to 
the unit sphere. Note / ij§ 2 G Cq'(S) and on we can determine 

i ? 2 (A §2 ) = fi^oY. 
Moreover, by Proposition 4.37, for r sufficiently small 

\\F 2 (f\ S 2)-f t :C 2 ' ri (S 6 ,h)\\<e/2. 

To finish the proof we make the following observations. Because of the uniform geometry of Sq in 
the h metric - for r sufficiently small- we can use standard linear theory on the interior to increase 
the L 2 norms of Appendix B to C 2,/3 norms. To make the dependence continuous, we let fi denote 
the normalized L 2 (S,h) projection of fi onto the span of ^(/jgz). 

The second half of the theorem is nearly identical. We sketch here a few of the main differences. 
Let S denote an S\p] and let Aj := A[p, ej,l], j = l,...,\E p \, denote the adjacent transition 
regions. Let Cj denote the center meridian circle of each Aj. As before S + denotes the connected 
component of S\ (u^lCj^j and AJ denote the other \E P \ components where the enumeration choice 
is the obvious one. Let Dj be the smallest geodesic disk in S 2 containing Y\p, ej, 1](A7). We set 

S := ({1} x 5+) U ({2} x D ± ) U ■ ■ ■ U {{\E P \ + 1} x D {Ep{ ) 

and g is the restriction of g§2 to the appropriate copy of S 2 . Let Y : S —> S such that Y_{x) := 
(l,y(x)) if x G S + and Y_(x) := (j + l,Y\p,ej,l](x)) if x G Aj . Again, using ip defined to be 
identically 1 on S and ip := ip[2d, d] ox_ on Aj for each j, we have the appropriate log cutoff function. 
From here, after defining F\,Fi in an analogous manner and compare with /j rather than fi, the 
result proceeds exactly as before. □ 

The extended substitute kernel. Following the general methodology of [14,18], we introduce 
the extended substitute kernel. The extended substitute kernel, which we denote by /C, is the direct 
sum of three-dimensional spaces corresponding to each standard region - central and non-central - 
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and each attachment. We first define the substitute kernel, the space of functions which will allow 
us to solve the semi-local linearized problem. 

Definition 5.20. We let 

IC v elCs:= K]p] K\p,e,n] 

peV(r) \p,e,n]eS{T) 

denote the substitute kernel. The spaces on the right hand side of the expression are defined below. 
Let /C[p, e, n] = (wi\p, e, n],W2\p, e, n], w$\p, e, n])]R where we define, for z = 1,2, 3, [p, e, n] G S(T), 
Wi\p + [e],e,n] := Ci\p+[e], e, n](^[2nP To(e) - 1, 2nP To(e) ](t))(^[2nP ro(e) + 1, 2nP To{e) ](t))f i , 

Wi]p-[e],e,n] := c i [p"[e],e,n](^[P[e] - (2nP To(e) - l),P[e] - 2nP ro(e) ](t)) 

■ty[P[e] - (2nP To(e) + l),P[e] - 2nP To{e) ](t))f i 

for each [p, e, n] G <S(r), where the coefficients Cj[p, e, n] depend smoothly on the parameters of the 
construction and are determined by requiring 

(5.21) / Wi[p, e, n]fj dh = 

J S\p,e,n] 

Notice the symmetry of the functions fj on S\p, e, n] in the h metric and the symmetric choice of the 
Wi\p,e,n] imply that the coefficients are not over deter mined. Moreover, the controlled geometry of 
S[p, e, n] in the h metric implies the coefficients are well controlled. 

To define IC[p] = {wi[p],W2[p],ws[p])r, we first need to determine a new cutoff function. Let 
Sym[p] C E>y p denote the largest region symmetric in all three functions x\,X2,xs such that 

dist§2(dSym\p], dS v ) > S. Because the graph is finite, for r small there is a uniform lower bound 
on the area of Sym[p\. In fact, under the assumptions of the construction, the distance between 
the centers of any two geodesic disks removed on M[p] must be at least tt/3 + S. The controlled 
geometry of Ey p implies that there exists a symmetric subset Sym'[p} C Sym[p] such that the 
exponential map from dSym'[p] into Sym'[p] is one to one on a 25 neighborhood. 

For X C S 2 , let T r (X) denote the r tubular neighborhood of X on S 2 where r is measured 
in the induced metric. Define the function ib\p] such that ip[p] = 1 on Sym'\p]\T2s(dSym'[p\) 
and smoothly cuts off to zero by dTs(dSym'\p\) n Sym'[p], remaining zero thereafter. Moreover, 
we can choose this cutoff so that it preserves the symmetry of the coordinate functions. That is 
Jg 2 i>\p\xiXj = if i ^ j. 

With this definition, we define the relevant functions on the central spheres. For i = 1,2,3, let 

Wi[p] := c-Y>[p]/j 

on S\p] where the coefficients depend smoothly on the parameters of the construction and are 
determined by requiring 

(5.22) / w i \p]fjdh = 5 ij . 

Js{p] 

Because of the construction, the coefficients are again determined solely when i = j as the other 
terms vanish; thus the problem is not overdeter mined. Moreover, the uniform lower bounds on area 
and the choice of domain where ip\p] is non-zero imply that the coefficients are uniformly bounded, 
independent of r. 

We now use elements of the substitute kernel to modify the inhomogeneous term and ensure that 
the modified term is orthogonal to the approximate kernel, see items (3), (4) below. 

Lemma 5.23. For each p G V(r), [p, e, n] G S(T), the following hold: 
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(1) Wi\p,e,n] is supported on S\p,e,n] and Wi\p] is supported on S[p\. 

(2) \\w t \p,e,n], Wi [p\:C 2 ^(M,x)\\<C. 

(3) For E G C°'P(S\p, e, n], x) there is a unique w G /C[p, e, n] such that -rj^E+w is L 2 (S\p, e, n], h)- 

orthogonal to the approximate kernel on S\p, e, n] . Moreover, if E is supported on Si \p, e,n], 
then 

\\w : C 2 ^(M,x)\\ < C(b,f3)\\E : C°^(S\p,e,n},x)\\- 

(4) For E G C°'P(S\p], x) there is a unique w G IC[p] such that ^jE+w is L 2 (S[p\, h) orthogonal 
to the approximate kernel on S\p\. Moreover, if E is supported on Si[p], then 

\\w : C 2 >?(M,x)\\ < C(b,P)\\E ■ C^(S[p], X )\\. 

Proof. The first and second items follow by definition of each of the functions and the uniform 
control on the coefficients, independent of r. The third and fourth items follow immediately from 
Proposition 5.17, the uniform equivalence of the h, x metrics on S±\p], S\ [p, e, n], and the definition 
of elements of IC\p], fC\p, e, n}. □ 

When solving the linearized equation, we must produce solutions on each extended standard 
region that satisfy sufficiently fast exponential decay conditions on one of the adjacent transition 
regions. (For central standard regions, the decay must hold on all adjacent transition regions.) 
This is the second use of the extended substitute kernel; we modify the solutions on C^[p, e, n"\ 
by prescribing the low harmonics of the solution to ensure such fast decay. The modification to 
the solution on each S[p], S\p, e, n] corresponds to an additional modification of the inhomogeneous 
term by elements of /C. To prescribe the low harmonics at each attachment to a central sphere, we 
introduce a three dimensional space of functions supported in a neighborhood of each attachment. 

Let IC\p, e] = (wi\p,e],W2\p,e],W3\p,e])M. where we define, for i = 1,2,3, \p, e] G ^4(T), 

(5.24) Wi ]p,e] :=c i [p,e}C h ^[0,l}oxV i [A[ P ,e,l},l,0}), 
where the Cj[p, e] are normalized constants so that 

cib,e]Vi[A[p,e,l],l,0] = l, 

(5.25) c 2 [ P ,e}V 2 [A\p,e,l},l,0} = cose, on C+[p, e, 0]. 

c 3 [p,e]V 3 [A[p,e, 1], 1,0] = sinfl 

The Ci\p,e] depend smoothly on the parameter r^(e). Moreover, Proposition 5.12 and (5.11) imply 
uniform upper and lower bounds on Cj[p, e] for i = 1, 2, 3. 
We now define the extended substitute kernel. 

Definition 5.26. Let 

JC := K v K S K A ■= K]p] K\p,e,n] JC\p,e] 

P eV(T) \p,e,n]eS{T) [ P ,e]eA(T) 

denote the extended substitute kernel. 

We introduce )C v \p, e, n], )C v \p, e] so that for v G K. v \p,e,n], Ch,v G JC\p,e,n] and likewise for 
tC v \p,e}. On each S\p, e,n], the substitute kernel itself can be used to modify the solutions to 
guarantee exponential decay. On the other hand, for the central standard regions we define the 
spaces /C[p, e] solely to arrange for the decay away from a central sphere. The prescription of the 
low harmonics near a central standard region is understood by considering the functions fi rather 
than fa. Recall that on C^\p, e, 0], the functions fi are multiples of v Ti ( e ) ■ e« where v Ti (e) refers to 
(3.7). 
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Definition 5.27. For \p,e,n"] G N+(T) and n" > 0, let 

]C v [p,e,n"] := (v 1 \p,e,n"],v 2 \p,e,n"],vs\p,e,n"]) R 

where Vi\p,e,n"] is supported on S\p, e,n"] and Vi\p,e,n"] := Ci\p, e, n"](fi\p, e, n"\ + m[p,e,n"]) 
where 

• fi\p, e, n"] are as in Proposition 5.17. 

• Wi\p,e,n"] is chosen such that — Chfi[p, e, n"\ + Wi\p, e, n"] is L 2 (S\p,e,n],h) orthogonal to 
the approximate kernel on S\p, e,n"]. 

• Ui\p,e,n"] solves the problem jChUi\p, e,n"] = —£hfi\p,e,n"] + Wi\p,e,n"] onS\p,e,n"] with 
zero Dirichlet data. 

• ci\p, e, n"] is such that on Cf [p, e, n"] we have ci[p, e, n"]f\ = 1, c 2 \p, e, n"]f2 = cos 9, c^\p, e, n" 
sin 9. 

Let fC v \p, e] := (vi\p, e], v 2 \p, e], vs\p, e])^ where Uj[p, e] is supported on A[p, e, 1] and 

Vi\p,e] := Ci\p, e]ip[0,l] o xVi[A\p, e, 1], 1, 0] 
where theci\p,e] are defined by (5.25). 

The choice of construction implies that a solution to the linearized problem on S\p, e, /(e)] will 
not be modified to induce the exponential decay along either adjacent transition region. On this 
extended standard region, the solutions will possess the linear decay guaranteed by Corollary 5.13. 
This will prove sufficient for the following reason. Solving the global linearized problem requires that 
the semi- local solutions coming from adjacent extended standard regions must have one solution 
decay exponentially on the common transition region (away from the nearest central sphere) and 
one decay linearly on the common transition region (toward the nearest central sphere) . The linear 
decay on both sides of S\p,e,l(e)] corresponds to the decay toward the nearest central sphere on 
each of its transition regions. 

Lemma 5.28. For £4 > there exists b large enough so that the Vi\p,e],Vi\p,e,n"],n" > defined 
previously are smooth on the corresponding S\p], S\p, e, n"] and satisfy 

(1) ChVi[p,e,n] G JC\p,e,n] and ChVi\p, e] G JC\p, e]. The support of e, n] is on S[p,e,n], 
and the support of C^Vilp, e] is on Si\p\. 

(2) Vi\p, e, n] = on dS[p, e, n], Vi\p, e] = on dS[p\. 

(3) \\vi\p,e,n] : C 2 ^(S\p,e,n},x)\\ < C(b,(3), \\ Vi \p,e] : C 2 ^(%], X )|| < C(b,(3). 

(4) For n" > 0, \\vi\p, e, n"] — 1, v 2 \p, e, n"] — cos 9, vs\p, e, n"] — sin 6* : C 2,/3 {Cf [p, e, n"], x)\\ < e 4 
and on C^[p, e, 0], v\ [p, e] = 1, v 2 \p, e] = cos 9, v^\p, e] = sin 9. 

Proof. Items (1), (2) follow immediately from the choice of construction. For [p, e], (3) and (4) 
follow from (5.25) and Proposition 5.12. For [p, e, n], we use the following. First, by Proposition 
5.17, \\Chfi[p,e,n] : L 2 (S[p, e, n], h)\ \ < e. Then the estimates from Lemma 5.23 for Wi[p,e,n] 
and standard theory imply ||uj[p,e,n] : L 2 (S\p, e, n], h)\\ < Ce. Upgrading the L? estimate to a 
C 2,13 estimate on the interior, and the uniform equivalence of the h, x metrics on S\\p, e, n], we get 
\\ui\p,e,n] : C 2 ^{Si [p, e, n], x)\ I < C(b,f3)e. As C x Vi\p,e,n] = on A ,i[p, e, n], A 1>0 \p, e, n + 1], we 
can apply Corollary 5.13 to get (3). Finally, (5.18) and the interior estimates on Uj[p, e, n] imply 
(4). ' ' ' □ 

Solving the linearized equation semi-locally. Using Lemmas 5.23, 5.28, we now solve the 
linearized problem on each extended standard region and obtain solutions with appropriate decay. 
Notice we presume the inhomogeneous term is supported in a small neighborhood of the standard 
region so the global proof must allow for precisely this setup. The proofs for central and non-central 
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(5.30) f[S](x) 



standard regions are essentially the same, though on each central standard region one must modify 
the solution on every adjacent transition region. 

Definition 5.29. Recalling Definition 5.19, on each extended standard region S[p, e,n] := S, we 
define the function f[S]: 

( i±i ,'frf A, 

l + ty(t), ifx = (t,0)eS 
e~ r ^, if x G A far ,n < 1(e) 

if x e Af ar , n = 1(e). 

When n < 1(e), take ip(t) to be a standard cutoff function that is 1 in a neighborhood of A c i ose n S 
and in a neighborhood of Af ar n 5. For n = /(e), take tp(t) = 1. 

Lemma 5.31. Let S 1 := S*[p,e,n] denote a standard region on M. Given 7 G (1,2), /3 G (0,1), 
there is a linear map 

TZg : {E G C°' /3 (S ; , x)-E is supported on Si} -»• C 2,/3 (S, x) x /C[p, e, n] 

suc/i i/ia£ t/ie following hold for E in the domain oflZg above and (p, w) = TZ^(E): 

(1) C x <p = E + l -0rw or C h p = ^E + w onS. 

(2) p vanishes on dS. 

(3) \\<p : C^(S,xJ[S})\\ < C(P,b,j)\\E : C°^(S 1 , x , f[S})\\ . 

(4) Fori = 1,2,3, \m\<C{b,P)\\E--C°' fi {SuX)\\ where w:=Y^ =1 ^iWi\p,e,n}. 

(5) TZg depends continuously on all of the parameters of the construction. 

Proof. Lemma 5.23 implies there exists w G /C[p, e, n] and ip' G C 2 ^(S) such that £ h <p' = ^E + 
w and <p'\gg = 0. The controlled geometry for r small implies (4) holds for where w := 
fJ-'iWi\p, e, n\. Now observe that the designated support of E and the construction of w implies 
that we can apply Corollary 5.13 to ip' on A c i ose to produce the necessary decay there. If n = 1(e), 
the same corollary provides necessary estimates on Af ar and the proof is complete. 

For n < 1(e), we consider how to modify <p' on Af ar to prescribe the fast exponential decay. 
Let ip' T denote the projection of ip 1 onto %i[Cf] where the of interest here is in Af ar . Let 
ip' ± = (p' — ip' T on C A f ar and let V^i := TZg(p' ± \ c +). Notice that (tp' — V <fi i)\ c + G %i[C^} and 
we denote 

(p' — V v >)\ c + = Qi + «2 cos 6 + 03 sin(9. 

Standard theory implies \ \<pt : C 2 ^ (S ,\)\\ <^C(b,P)\\E : C°' p (S,x)\\- Coupling this with Corollary 
5.9 (4) implies \o>i\ < C(b,f3)\\E : C 0,/3 (S,x)\\- We now choose v G K. v \p,e,n] uniquely such 
that on C± , ip' + v = + T^d(v±\ c +). We may write v = Y^,=il JL i v i\P-> e i n ]- For each i, let 
V{ := 1Zd(vi\p,e,n]±\ c +). It is enough to show that for each i = 1,2,3, there exists a unique /// 
such that 

A^Yi ~ v i)\c+ = a i( S n + cos65 i2 + sin0<5j 3 ). 

Lemma 5.28 (4) and Corollary 5.9 (3) give uniform (positive) lower bounds on the difference of the 
coefficients of Vi — Vi at Cf. This, coupled with the bounds on the \ai\ implies 

\tf\<C(b,mE:C°*{SuxJ[S\)\\- 

Setting ip = ip' + v and w = w + ChV, we see Ttgp^_ = ip on and thus Corollary 5.9 provides the 
necessary decay on Af ar . □ 
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The procedure for solving the linearized problem with appropriate decay at each central standard 

region is nearly identical. In this case, we use IC\p, e] rather than /C[p, e, n] to modify our solutions 

but all of the estimates will follow by the logic outlined in the previous proof. 

\A\ 2 

Observe that in the statement of the lemma, the term h^r does not appear. This is a simple 
consequence of the fact that p 2 = \A\ 2 /2 on supp (w) C M[p]. 

Lemma 5.32. Given 7 £ (1, 2), /3 G (0, 1), there is a linear map 

lZ~ s[p] : {E G C°>P(S\p],x) : E is supported on S^p}} -> C 2 ^(S[p], X ) x (/C[p] U eeEp K]p,e\) 

such that the following hold for E in the domain ofTZg^ and (<p,w) = IZg^(E): 

(1) C x ip = E + w or Ch(f = ^ApE + w on S\p]. 

(2) ip vanishes on dS[p\. 

(3) \\<p : C^(S\p], X )\\ < C(b,P)\\E : C^(Si[p],x)\l 

(4) EMM\ + ^ij\^\P^ j }\<C(b,P)\\E-.C^(S 1 , X )\\ where 
w = Y,i Vi\p\ w i\p\ + Eij A*i[P. ejH[p, ej]. 

f5j ||<p : C 2 'P(A[p,e, l],x,e"^)|| < C(b,(3, 7 )\\E : C^(S 1 \p], x )\\ for all e G £7 p . 
ftfj ^g[p] depends continuously on all of the parameters of the construction. 

Proof. Lemma 5.23 implies there exists w G /C[p] and G C 2 ' /3 (S'[p]) such that -C/^' = pp-E 1 + w 
and v'lggjp] = 0- For each A[p, ej,l] where ej G -E p we modify <p' in a manner identical to the 
previous proof. On each Cf\p, ej, 0], find the unique Vj G /C„[p, ej] such that for ^ = ip' + ty, and 
¥>j,± representing the part orthogonal to T-L\[Cf\p, ej, 0]] on C^[p, ej,0], one has 1Zg((pj : ±) = ipj. 
The estimates on tpj , i>j and the unique choice of Vj follow as they did previously. Now set 

\Ep\ 

3=1 

Observe that supp {v{) n supp (vj) = for i ^ j and thus f\A[p,e-,i] = fj- Finally, setting w = 
w + Y2j=i £hVj gives the result. □ 

Solving the linearized equation globally. We can now solve the linearized problem semi-locally 
on each transition region and each extended standard region. To get the estimates we desire, the 
inhomogeneous term for each extended standard region must be supported in a neighborhood of the 
standard region. Thus as a first step toward solving the global linearized problem, we determine a 
partition of unity on M that we apply to the inhomogeneous term. Solving the problem on each 
region separately, we then patch the solutions back together to construct a global function. Of 
course, this process will introduce some error, which we correct by iteration. 
We first introduce all of the cutoff functions we require. 

Definition 5.33. We define uniquely smooth functions ^s\p]^s\p]^S\p,e,n],^s\p,e,n]^^\p,e,n'] such 
that 

(i) i>s\p] = i>s\p] = 1 on s \p\- 

(ii) ips\p] = V^l) 0] o x on each A\p, e, 1] for e G E p . 
(Hi) V^jp] = V^O; l]°x on each A[p, e, 1] for e G E p . 

M 1pS]p,e,n] = ^S\p,e,n\ = 1 0H %» e > n ]- 

(v) For n < 1(e), ip S \p,e,n] = ^[1,0] o x on A[p,e,n] and ^s\p,e,n] = ^[1,0] o x on A\p,e,n + 1]. 

(vi) ips\p,e,l(e)] =ip[l,0]ox on A[p,e,/(e)] and ip S \p,e,l(e)] = Y>[l,0]ox on A[p~[e],e,l{e)\. 
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(vii) For n < 1(e), 4>§[p,e,n] = ^[°> l ]°xon A[p,e,n] and i/>s\p,e,n] = V'fO, 1] ° x on A\p,e,n + 1]. 
(OTMj ^^^(e)] =#,l]°i on A[p,e,Z(e)] and V5[ P ,e,/( e )] =#,!]« on A[p"[e], e, Z(e)]. 
^ ^A[p,e,n'] = 1] ° I on A[p, e, n']. 

Observe that ifrsWAi V\S[p,e,n]> ^A[p,e,n'l form a partition of unity on M. Also note that each of 
the functions ^sM'^Slpen] are identically 1 on almost all of S\p], S\p, e, n], respectively. Near the 
boundary they transition smoothly to zero. Finally, supp(tps[ p ]) C <Si[p], supp(ips[p,e,n]) C <Si[p,e,n]. 

We now set the notation for defining a global C 2 ' 13 function by pasting together appropriately 
cutoff local functions. 

Definition 5.34. Let u[p] G C k ^(S\p}), u\p,e,n] ^C k ^(S[p,e,n}), p G V(T),\p,e,n] G S(r), 6e 
functions that are zero in a neighborhood ofdS[p], dS[p, e, n]. We define U = \J({u[p], u\p, e, n]}) G 
C k,/3 (M) to be the unique function such that 

(i) U\ S [p] = u\p],U\s\p, e ,n] = u[p,e,n]. 
(H) U \A\p,e,l] =u[p] + J2eeE p u \P> e i 1 ]- 
(Hi) Forn' < 1(e), U\ A ^ e>n >] = u[p,e,n' - 1] + u[p,e,n'}. 

(iv) ForU\ K[p+[e]te ^ e)] = u[p + [e],e,l(e)-l]+u[p + [e},e,l(e)] while £/|A[p-[ e ],e,;(e)] = u \p~ N> e > K e ) ~ 
l]+u\p+[e],e,l(e)}. 

We define the global norms that we will use in the statement and proof of the global linearized 
problem and in the main theorem. 

Definition 5.35. For k G N, (3 G (0,1), and 7 G (1,2), we define a norm \\ ■ \\k,p n on C k ^(M) 
by taking \ \f\\k,p,^ to be the supremum of the following semi-norms with p G V(r), [p, e, n] G S(T), 
and \p,e,n'] G N(T) defined in 4-8: 

(i) H/rC^SMx)!! 

(ii) r-^||/:C fc ^(S'i[ P ,e,n],x)|| 

(Hi) T ~^'^)\\f : C k ^(A\p,e,n'],x,e-^)\\ 
For w G KZ such that 

W = ^Z\ ^2 wbl'filP] + Vi[p,e]wi\p,e}+ ^ /J-i\p,e,n\wi\p,e,n] J , 

»=i \pev(r) [p, e ]eA(r) [p,e,n]eS(r) / 

Zei nj x denote the projection of w onto KZx ■ We define norms for each of these projections separately 
to help highlight the differences in scale for the different coefficients. To that end, let 

(5.36) := max ||u> A ||a := max {\m\p, e]|}, ||«; s || 7 := sup {T~ in \m\p, e,n]\}. 

1=1,2,3 «=1,2,3 i=l,2,3 

p [P,e] [p,e,n] 

We are finally ready to state and prove the main result of this section. The strategy of proof 
is to first use Proposition 5.8 to solve the problem on each transition region. We then apply the 
cutoff functions to the inhomogeneous term and use Lemmas 5.31, 5.32 to solve semi-locally on 
each extended standard region. Finally, we patch together these semi-local solutions via 5.34. This 
process introduces error which we show is small enough to correct by iteration. 

Proposition 5.37. There is a linear map Km ■ C°>"(M) -»• C 2 'P(M) x K such that for E G 
C°>P(M) and (tp,w) = TZ M (E) the following hold: 

(1) £ x ip = E+^w on M. 

(2) |M| 2 ,/3, 7 + \\w v \\v + \\w A \\ A + lk 5 || 7 < C(b, /3,7)l|£||o,/3, 7 

(3) TZm depends continuously on all of the parameters of the construction. 
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Proof. We first prove the proposition, assuming supp (E) C (Uy(r)Si[p]) U (l)s(T)Si\p,e,n]). Using 
Lemmas 5.32 and 5.31 we determine pairs 

{(p\p,e,n],w\p,e,n\) = ^^(^Is^n]), = n s\p]( E \si\p]) 

for each \p, e, n] G S(T),p G V(r). (Note that w\p] will be a linear combination of elements of both 
IC[p] and JC\p, e].) In order to patch together the local solutions </?[p], <p[p, e, n] and maintain the 
regularity, we have to cutoff each of these functions near the boundary. But this introduces some 
error. We denote 

TIE := \J({^ [p] ^\p],^ [pte ^\p,e,n]}) G C 2 ^(M), 

\v(r) S(r) / 

EE := U({[^ [p] ,C x Mp], {^~ s[p ^ nV C x MP^M}) e C°^(M), 

where here [ , ] denotes the commutator. That is, [Y^p Ad^tp] = ^sip^x^lp] ~ ^xW'sjpjV'tp]) an( ^ 
the like for <p\p, e, n\. 

One easily checks that, as the support of E implies i/}g^jZ x ip\p] = C x <p[p] and the like for 

ip\p,e,n], 

\A\ 2 

(5.38) C X TZE + 8E = E + i^WE on M. 

Moreover, for 7' G (7, 2), we get the estimates 

\\1lE\\ 2Al <C{b,[} n )\\E\\ ^ 
\\WE v \\ v + \\WE A \\ A + \\WE s \\ y <C(b,p,j)\\E\\ 0Aj 

wsEWo^Kr^-^c^p^^'mwo^. 

The first two estimates are immediate from Lemmas 5.32, 5.31 and the uniform geometry of the 
immersions on the support of WE. To determine the third, we first observe that applying Lemma 
5.31 to EE for 7' G (7, 2) yields the estimate 

||^||o,/3,7<e {7 - 7 ' k C(6,/3,7)||ii;|| OA7 . 

Next observe supp (EE) C (u peV ( T Si\p}\S\p\) [j ( ^[p,e,n]eS(r)Si[p, e, n]\S[p, e, nfj . Finally, for a 
fixed e G E(T) U R(T), on the supported regions, by (5.4), we have the estimate 

exp(( 7 -7'k) < C(7,7',&)|rd(e)| (y - 7) < C( 7 , 7, &)t ( ^. 

If r is sufficiently small, we can guarantee C( 7 , 7', fyr^'^^ < 1/2. This allows us to remove the 
error by iteration. 

Now observe that since EE G C ' 13 (Uv(r)Si[p]) U ( u 5(r) &i \p, e i n \)-, we can apply the steps above 
to EE to produce 7l}E, W 1 E,£ 1 E. We then proceed by induction to produce, for all r G 

\A\2 

(5.39) £ x ^ r ^ + £r+lE = £TE + ^2 W ' E on M - 
Defining 

00 00 

(5.40) <p := K r E; w := W r £, 

r=0 r=0 

which converge based on the provided estimates, we prove the proposition in the first case. 

40 



Now consider the more general case. We begin by applying Proposition 5.8 to each A[p, e, n'\ to 
determine 

V[p, e, n] = fcA\p, e ,n']E\A\p,e,n'] 

such that V\p,e,n']\c+\p je ,n'-i] £ / Hi[C + \p,e,n' — 1]] and V\p, e,n']\c-\p, e ,ri] = 0- We define 

E := \5{{^ sb>] E^ s ^ eM E}) + U({0, ^a^,], £ x ]F[p, e, n']}) G C°^(M), 
and observe that i£ C (Uy(r)£i[p]) U i}^S{T)S\\p, e, n]). We now apply this case to the work done 

2^ 



above to produce (<p, w) G C 2,/3 (M) x /C such that £ x </3 = on M. We define </9 G C 2 ^(M) 



such that 

= ip + V({0,^ A[Ptey] V[p,e,n}}). 

Recalling that ips]p] > V'S^e.n] > "0A[p,e,n'] partition M, we see £ x <£ = + ^-ro. Moreover, 
satisfy the necessary estimates so the proof is complete. □ 

6. The Geometric Principle and proof of the Main Theorem 

Throughout this section, let 7 G (1,2),/? G (0,1), 7' G (7,2), and /?' G (0,/3). Any constant 
depending only on b, j3, 7, /?', 7' we will simply denote by C. If it depends on other constants, they 
will be made explicit in the notation. 

Prescribing the extended substitute kernel on a non-central standard region. As in 

the paper [6], prescribing the extended substitute kernel on the non-central standard regions can 
be done entirely on the linear level. The main tool of the argument is Green's second identity 
[4], which is a linearization of the force balancing. Before we proceed, recall that fa, defined in 
(5.15), are simply translation vector fields that have been projected onto the normal bundle of 
(multiplied by some constant scale factor). Since translations preserve the mean curvature, Cfa = 
for i = 1,2, 3. 

We begin by defining, for i = 1,2, 3, 4>\\p, e, n] G C 2,/3 (S[p, e, n\) where 

<£i[p,e,n] = V[0, 1] °xVi[h c i ose ,ci[p,e,n]l Kclosa ,ti\, 
(6.1) <j)' 2 \p,e,n) = ip[0, 1] oxV2[A c i ose ,c 2 \p,e,n\sinh£ Aclose ,0], 

^b» e . n ] = ^[0,1] oxV 3 [A dose ,C3\p,e,n\ sinh &\ close , 0] . 
We choose Cj[p, e,n], i = 1, 2, 3, such that - recalling (5.21)- 

(0.2) / M|p]£._|> Ae ,„ ]d9 = _ % 

^ C 'dose 

where here C c \ ose represents the middle meridian circle on A c ; ose . (Here A c / ose ,Aj ar are as in 
Definition 5.19.) For r d (e) > 0, (3. 6), (3. 7) imply that dfj/dt = on C dose for j = 1,2,3. For 
Trf(e) < 0, the same equations imply that f 2 , fo = 0, dfi/dt = on this meridian. 

For the rest of the argument we drop the dependence on [p, e, n] G S(T) and presume the reader 
understands the implied dependence. Where there might be confusion, we state it explicitly. 

Notice that supp(<^) C A c / ose and supp (£ x <^) C SinA c i ose . Proposition 5.12, (5.11), (5.4), and 
the definition of fa imply 

C 1 - 1 < Ci < C, i = 1,2,3 

and 

||# : C 2 ^(A close ,x,e-^\\ < Ce e < Ct~ 1 
||£ X ^:C^(S, X )||<C< 
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Observing supp {C x <j>'j) C S\ PI A c ; ose , let ((f)", = 7£g(— £ x <//) and finally, define 

0i := <Pi + <t>'i- 

Lemma 6.4. Given 4>%,w_i as defined above, for S := S\p, e,n], 

(1) C x 4>i = l -0rw { on S. 

(2) \\^:C 2 ^(S\A close , X )\\ <C. 

(3) \\<j>i : C 2 >P(A close ,x,e-^)\\ <Ct~\ 

(4) \\(f>i : C 2 ^(A /ar , X ,e-^)|| < C when n < 1(e). 

(5) Ufa: C 2 >P(A fariX , T^)|| < C when n = 1(e). 

(6) \w i — Wi\ < C/\ logr| where Wi represents the appropriate Wi\p,e,n]. 

(7) 4>i,w.i are unique by their construction and depend continuously upon all of the parameters 
of the construction. 

Proof. Items (1) and (7) follow immediately from the definition. Items (2) and (3) follow from 
(6.3), (5.4) and (4), (5) follow from the definition and Lemma 5.31. 

For (6), we require the use of Green's Theorem. For a domain Q C S, we have 

L /i£ ** - ^ Ji dx = L h (W>- * (f ) ' 

where here r\ is the outward pointing conormal to dQ. Let be the connected component of 
S\ (c dose U Cf ar ^j where here Cf ar is the middle meridian circle on Af ar . As translations preserve 

\A\ 2 

the mean curvature equation, we know C x f\ = 0. Further, using C x 4>i = Li-Wj, the above equation 
reduces to 

L i* ^ - L - - L 0# - $*) ~ L - W • 

Using (6.2), for r^(e) > 0, since dfi/dt vanishes on each boundary curve it is enough to show 
appropriate control on Recalling the construction of (f>", observe that C x <t)'- = on Af ar and 
A c / ose \ (S 1 n A dose ). If we denote 0S l n A c/ose := C~ ose l then ||^' : C 2 ' /3 (C+ r U Cj Me>1 , M)\\ < C. 
Item (2) in Corollary 5.13 then provides the necessary estimate. 

For Td(e) < 0, the same argument works for f\. For f2,fs, we observe that on C c i ose ,Cf ar , 
\9f2/dt\, \dfs/dt\ < C\Td(e)\ l l 2 . Since fi = = on the boundary curves, using the estimate for 
<j>" on C^ ar U C~ lose 1 coupled now with (1) in Corollary 5.13 is strong enough to provide (6). □ 

Prescribing substitute kernel on the central spheres. The prescription of the extended 
substitute kernel on each central sphere is completed in two steps. We first prescribe substitute 
kernel based on the unbalancing parameter d. We then prescribe the extended substitute kernel 
at each attachment based on the dislocation parameter £. Part of the difficulty in what follows 
comes from the fact that the unbalancing is understood by considering the graph F(d, 0) and an 
immersion based on that structure while the dislocation is understood by considering dislocations 
from the graph T(d, £), as the dislocations are introduced on an immersion positioned relative to 
that graph. 

We state here an important lemma without proof. The lemma depends upon a flux argument 
[21, 22]; the proof for our specific setup can be found in Appendix A. 3 of [10]. We state the lemma 
in the form exactly as we need for our purposes. Throughout this subsection, for any p G V(r), let 

S + [p] denote the connected component of S[p]\ (^J ee E p C[p, e, 1]^ . 
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Lemma 6.5. For d,£, Y d £ and the function H g i u ing induced from this immersion, 
/ H g i uing \p] N dX dg = 7rY] T d (e)0- e ,pR'[e]ei := Trd c (p). 

Here p' £ V(T(d, £)) is the vertex corresponding to p and e' corresponds to e. 

Notice we do not expect that d equals d^. However, the difference between the two is bounded 
by higher powers of r. 

Proposition 6.6. Given d, £, let H g i u i ng denote the gluing function induced from the immersion 
Y d:C . For eachpe V(T), there exist (\) g i uing \p\ G C 2 ^(S\p\), £t- [p] , //■ [p, ej], i = 1, 2, 3, j = 1, . . . , \E p \ 
such that 

( V t-X^gluinglp] + P~ 2 H g luin g = E rfM W i\p] + E Kb' e j\ w AV, &j\ 0n ^\p\ ' 
(2) 4> g luing[p] =0 on S[p\. 



\^ g imn g [p}:C 2 ^(S[p],x)\\<Cr. 

\<P g i wng [p] : C 2 ^(A[p,e,l],x,e-^)\\ < Ct for all e G E p . 

E?=iMj[p]ei-d(p)| <er/16. 
KM|<Cr. 

f7j E'ac/j of the (j> g i U ing, Kb]; Kb> e j] are a ^ unique by their construction and depend continu- 
ously on all of the parameters of the construction. 

Proof. Determine K sucn that for j = 1, 2, 3, 



(3) 
(4) 
(5) 
(6) 



L [ J2 Vi Wi ^ - TJ\2 H gluing I /j [p] = 0. 



Rather than investigating this integral, we initially consider 

3 



^ ^ ^X^K^b] - |^|2^ 9 /mn 9 ^ N d< dh, 



recalling N d ^ • = 7r/j. Based on the definition of iOj[p], 



An application of Lemma 6.5 implies 



where d^(p) is defined in Lemma 6.5. By Proposition 5.17, \\fi\p\ — fi\\c k (S 5 [p],h) ^ e ) where e can 
be made as small as we like by making r small. Coupled with Proposition 4.46 and Lemma 5.23, 
we estimate 

3 , . 9 \ , 

dh 



J s+[ ] [ibtivifc] - J^2 H 9luingj (/>] " fj) 



< Cer < er/32. 



Notice C is independent of r so we can easily achieve the inequality Ce < e/32. Observe also that 
d c (p) - d T{d/) {p) + d r(d/) (p) - d{p) = ^2 T d( e )o- e , P (R'[e]ei - R'[e]ei[e] + R'[e]ei[e] - vi[e;d,0]) . 
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Using (4.21) and Proposition 4.23 we see 

\d c (p)-d(p)\ < CtV 2 + CCt 2 

and item (5) follows. 

The rest of the estimates follow immediately from Lemma 5.32 and the bound on p~ 2 H g i u i ng . □ 

Prescribing the extended substitute kernel on the central spheres. We are now ready 
to describe the prescription of elements of /C[p, e] on each extended central standard region. The 
proof relies on the estimates of Appendix A and uses the functions defined in Appendix B. Given 
the parameters d, £, we modify the inhomogeneous term of the linearized equation by elements 
of ®ee_E p ^[Pi e ] such that the coefficients of these elements are close to the parameters C[p> e ]- 
Notice that when we solve for p~ 2 H dislocation,! the elements of the substitute kernel have coeffi- 
cients bounded by C||£||/| logr| which reflects the fact that C[p> e ] is essentially orthogonal to the 
approximate kernel. 

Definition 6.7. We normalize the functions fi defined in (5.16) on the meridian circle Cf [p, e, 0]. 
Precisely, for [p,e] G ^4(r), choose d { \p, e] such that on Cf[p,e,0], 

ci[p,e]7r/i = 1; c' 2 [p, e]7r/ 2 = cos6»; d 3 \p,e]Trf 3 = sinfl. 

Let 

(6.8) c := max |c'ta,el|. 

v ' i=l,2,3, 

MeA(r) 

The normalization implies that on Cj^fp, e,0] the function that describes as a graph over 
Y + C [p, e] is approximately Ci [p, e] / d 1 [p, e] + (( 2 [p, e] /c' 2 [p, e] ) cos 9 + (£ 3 [p, e] / c' 3 [p, e] ) sin 9. (See 
Appendix B.) Observe that as r — > 0, c' converges to a value that depend upon the geometric limit 
(in this case § 2 ). That is, d depends on b and is independent of r. 

Proposition 6.9. For each p G V(r), i/iere exists 4> dislocation^} G C 2 '^(<5b]) anc ^ a4'[p]> a4'[p> e j] 
with i = 1,2,3 and j = l,...,|£p| such that, for parameters d,C, the immersion Y^, and the 
induced function H dislocation, 

(1) C x (f)di s locaUon[p\ + P' 2 H dislocation = Ya=1 (^"[pHM + Y^fll PilPi e j] w i\Pi e j]j on S[p\- 

(2) <\> dislocation [p] = 071 SS'fp]. 

W + ICiM/cJfoe,-] - Pi\P^j\\ < C\\C\\/\ logr| /or i = 1,2,3 and each ej G £„. 

(4) 110 dislocation [p]:C 2 ^(S[p],x)\\<C\\C\\. 

(5) 1 10 dislocation 

[p] : C 2 ^(A[p, ei ,l], X , e -^)|| < C||C||/|logr| /or eac/i ej G £7 p . 

^dislocation 

[p], /j'/[p], /j'/[p, ej] are all unique by their construction and depend continuously on 
the parameters of the construction. 

Proof. Throughout the proof, i G {1,2, 3}, j G {1, . . . , \E p \}. Let f[p] represent the function 
from Definition B.5. We construct dislocation^ G C 2,l3 (S[p\) in the following way. On M[p], 
fet <t> d i s i oca ti n\p} = f\p]- 0n M H n %]> for P = P + M set & : = ^% + 3, a + 2](t) and for 
p = p~[e,] set fj := i/>\P[ ej ] - (a + 3),Pfo] - (a + 2)](t). On M[ej] D %], let 0^ ocation [p] = 
$j'/[p] — (1 — V , j) 7r Si=i Ci[P) e j]/i- Finally, on each A[p, e.,-,1], find such that C x V_j = 0, 
Zj = -7T £- =1 C»[p, ej]/i on C~ [p, ej, 1] and = on C+[p, e^, 0]. On A[p, ej , 1], let <j)' dislocation [p] = 
— 7T X)i=i Ci[p> e j]/« + (1 — ^S[p])3^j- The construction of implies 

HZj : C 2 ' /3 ([o,6 + 2] x S 1 n M[ej],x)|| < C||CII/I lo g r l- 
Noting the estimates provided by Corollary B.6, we have the following estimates for $ dislocation \p\ ' 
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(!) ^'dislocation^] is supported on S x \p] and fi dislocation \p] = on 

(2) For each j, \\C x< f/ dislocation \p\ - p~ 2 H dtslocaUon : C°^(5i[p] n M[e,], x )|| <C||C|/|lo g 



T . 



Recall the argument used in Lemma 5.31 to determine the coefficients /Xj. In this case we 
determine /^[p, e j] by noting that on [p, ej, 0] we must have 

This is a simpler expression than the one in Lemma 5.31 because Vi,f\ G %i[C^]. Then 
Corollary 5.9, items (3) and (4) coupled with the estimate on V_j imply 

ICi[p,ej]/cJ[p,ej] -fi.\p,ej]\ < C\\C\\/\ logr|, and 

3 

\Wdisiocation\p\ + J>>' e >M ■ ^(Afc, e,-, 1] , x , e"^) 1 1 < C||C||/|logr|. 

i=l 

Using Lemma 5.32 with E := ^'dislocation^] - p~ 2 H dislocation, let 

{ftdislocationlpli w dislocation) ■= ^5[ p ](- E ') 

where 

Wdislocation = ^ Vi\p\ w i\p\ + e A W i\P' 

i i,j 

Then: 

( 3 ) ^"dislocation^] = ^'dislocation^] ~ P~ 2 H dislocation + Wdislocation On S[p] 

(4) |<W|,K"|p,e i ]|<C||C||/|logr|. 
Set 

4>dislocation\p] = <P'dislocation\p] ~ 4> dislocation^] ~ Yl Hi b> e # b> e i] 

i 

and 

/4'[p, ej] = n"'\p, ej} - (ijp, ej}. 
We complete the proof by appealing to all of the estimates above and those of Lemma 5.32. □ 



Prescribing the extended substitute kernel globally. Let M ,M denote the finite dimen- 
sional vector spaces with components on M and indexed over i = 1, 2, 3 and p G V(T) or i = 1, 2, 3 
and [p, e] G -A(r). Let M s denote the infinite dimensional vector space indexed over i = 1, 2, 3 and 
[p, e, n] G S(T). Let £ G M y x R A x M 5 and let £ x denote its projection onto R x . We define norms 
such that 

H^llv :=max{|&[p]|}, \\^ A \\ A := max{|&[p, e]|}, ||£ 5 || 7 := sup {r^^p, e, n]|}. 

l 'P 1 >M i,[p,e,n] 

We proceed by assuming throughout 

(6.10) \\Z v \\v <st/V3, \\£ a \\a<Ct/(/, ||£ 5 || 7 <t. 
Now for p G F(r) and [p, e] G A(T), set 

(6.11) d(p) = ^&[p]e; and Ci[P> e ] = c 'iiP, e ]&[P> e ]- 

i 

Notice the scaling for the norms for the spaces M y ,IR A implies that ||c2||d < £T, \ \C\\ < Qr. With 
this definition of d,C we find <\> g i u ing [p] j ^dislocation [p] from Propositions 6.6, 6.9 and <^[p, e,n] from 
Lemma 6.4 and set 

^ = U ^5 b] • i<Pgluing[p] + <P dislocation 

i I 
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Setting fj,i\p] := /i/'Jp] + p![ [p] , m [p, e] := /u'Jp, e] + fjf-\p, e] where these coefficients come from Propo- 
sitions 6.6, 6.9, define 

— := ^2 f J 'i\p\ w i\p]+ m\p,e]wi\p,e]+ ^2 &[p,e,nfe[p,e,n] € K. 

i=i \pev(r) MeA(r) [p, e ,n]eS(r) / 

Here Wj[p, e,n] comes from Lemma 6.4. Using Proposition 5.37, determine 
Now set 

$1 := + and ■= w" + u/. 

Proposition 6.12. anc? as defined above for the immersion where d,C are defined by 
the components chosen in (6.11) ; depend continuously on £ and satisfy: 

(1) £ x $ € + p- 2 H error = on M. 

(2) ||^|| 2i/3i7 <C(||^|| 7 + r+||C||). 

(3) \\{w $ -w $ ) v \\ v <er/8, 
\\{w^-w^) a \\ a <Ct/c\ 
||(^-^) 5 || 7 <C||^|| 7 /|logr|. Here 

w i'=^2i\ ^2 &\p\ w i\p\ + €i\p, e ] w i\P> e ] + ^ &[p, e, n]wj[p, e, n] J . 

i=l \peV(r) [p,e]6A(r) [p,e,n]eS(r) / 

Remark 6.13. Notice that the definitions of the norms implies H^IU = H^Ha, ||£^||v = H^I'llv') 
and | \£ s \ | 7 = | liy^l | 7 . 



Proof. From the definitions and previously determined estimates, we observe C x ^£. + p H error = 
\0-W£ on M and ||^|| 2 ,/3, 7 < C(\ \£ s \ | 7 + r + ||C||). Moreover, as the function -C x ^ + - 

-2] 



p 2 # error is supported on U peV ( r )Si [p^-S^] (J U\p : e,n]eS(r)Si b> e, n]\5[p, e, n], we can apply the 
same technique as in the proof of Proposition 5.37 to determine 

|| - £ x &t + ^w' - p- 2 H error \\ 0Al < Cr^-T(||^|| 7 + r + ||C||). 

Finally, observe that a term by term comparison implies 

\\(u/ -w^) v \\ v < er/12 
\\(w-w^ a \\ a <Ct/c' 

\\(v/-we)%<c\\z%/\]ogT\. 

□ 



7. The Main Theorem 

We are now ready to state and prove the main theorem of the paper. The proof of the theo- 
rem depends upon all of the estimates in the paper, but the main estimates we need come from 
Propositions 6.12 and 5.37 as well as the estimates in Appendices A,B. 

46 



Theorem 7.1. Let T be a flexible, central graph. There exists C > sufficiently large and r' > 
sufficiently small such that for all < r < r': 

T/iere exist C G K A and d G £>(r) mtt ||C|| < CV, ||d|| D < er, and / G C 2 '^(M) such that 
{Yd,c)f '■ M — > M 3 is a well-defined immersed surface of constant mean curvature equal to 1 and 
ll/l|2,/3,7 < r 15 / 16 . Moreover ifT is pre-embedded then (ld,c)/ * s embedded. Here M is the abstract 
surface defined in (4.7) depending on V and r, * s ^ e initial immersion described in Section 4 
depending on T and the parameters r,d,C, and (Y^)/ is defined in Appendix A. 

Proof. Choose C so that C > 2C for all C appearing in the statement of Proposition 6.12. Choose 
t' small enough such that the value max e6E / r \ uR ( r ) |r'f(e)| satisfies all geometric requirements 
throughout the paper. Notice that none of the geometric estimates t' must satisfy depend upon 
the structure of a graph but only on the function f . For any < r < r', let 

B :={u G C 2, P(M) : |M| 2 ,/3, 7 < r 3/2 }x 

{£eR v xR A xR s : \\£ v \\ v < er/V3, \\£ A \\ A < Ct/c', ||£ 5 || 7 < r}. 

We define a map J : B — > B in the following manner. Let (n, £) G B. Recall £ v induces a 
d G D(T) and £ A induces a C G M A as outlined in (6.11). With this d, C, we find a F(d,£) G 7"(r) 
and determine the immersion Y d ^. Set v := — «. Then by Proposition 6.12, C x v = L-L-t^ — 

p' 2 H error - C x u and | |u| | 2j/3 , 7 < C(| | 7 + r + | |C| |) + r 3 / 2 < r 15 / 16 . Now apply Proposition 5.37 
to obtain (u 1 , w') := 1Zm{p~ 2 {H v — H^) — C x v) where here H v is the mean curvature of the surface 
induced by (Yd^) v . Then 

(1) C x u' = p- 2 (H v - H dtC + H error ) + + C x u 
and by Proposition A. 5, 

(2) \\u'\\ 2Al + \\w% < Cr 15 /8 r -V8 < r 3/2. 

Define /i G R v x M A x such that 

^2 fJ-i\p\ w i\p\ + ^2tH\p, e]uii[p, e] + ^2ni\p,e,n]wi\p,e,n] = w' + - w^. 

By the estimates of Proposition 6.12 and the equivalence of the norms of fi and w^, 

(3) \\n V \\v < er/V3, \\n A \\ A < Qt, ||/x s || 7 < r. 

Define the map J(u,£) = {u',n). Then by (2), (3) J{u,£) G B and the map J : B ->• £ is well 
defined. Moreover, for some ft G (0, /3), B is a compact, convex subset of C 2 '^' (M ) x (M y x R A x M s ) 
and one can easily check that J is continuous in the induced topology. Thus, Schauder's fixed point 
theorem, [4] Theorem 11.1, implies there exists a fixed point (u' ',/*') G 5. Let / = 3> M / — n'. 

Since (u', n') is a fixed point, the definition of /z' implies w' — = 0. Moreover, C x u = C x v! . 
Thus 

Hf — Hd^ + H error = 0, or Hf = — H error = 1 
and the surface (Yd^)f has mean curvature identically 1. Embeddedness follows when T is pre- 
embedded as in this case M dX is embedded and ||/||2,/3, 7 < r 15 / 16 . □ 

Appendix A. Quadratic estimates 

For completeness, we include here a proposition we will need, the proof of which appears in the 
appendices of [6,20]. Let X : D — > U be an immersion of the unit disk in M 2 into an open cube 
U CM 3 equipped with a metric g. Assume dist g (X(D), dU) > 1 and there exists c\ > such that: 

(A.l) \\8X :C 2 ^(D,g )\\< Cl , \\ gij ,g^ : C 4 ^(U,g )\\ < a, go< Cl X*g, 

where here dX represents the partial derivatives of the coordinates of X, g tJ are the components of 
the inverse of the metric g, and go denotes the standard Euclidean metric on D or U respectively. 
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We note that (A.l) can be arranged by an appropriate magnification of the target, which we will 
exploit in order to make use of the following proposition. 

Let v : D — > R 3 be the unit normal for the immersion X in the g metric. Given a function 
4> : D — > R which is sufficiently small, we define : D — > U by 

(A.2) X^p) := exp x(p) (0(p)i/(p)) 

where here exp is the exponential map with respect to the g metric. Then the following holds: 

Proposition A. 3. There exists a constant e(c±) > such that if X is an immersion satisfying 
(A.l) and the function <f> : D — > R satisfies 

||^:C 2 ^(A5o)||<e( Cl ) 

then Xff, : D — > U is a well defined immersion by (A.2) and satisfies 

WXj, — X — (f)u : C^{D,go)\\ < C( Cl )||0 : C 2 >? (D , g Q )\\ 2 

and 

WH4 -H- C x * g ■ C°>P(D,g )\\ < C{ Cl )\\<f> : C 2 >?(D, g )\ | 2 . 

Here H = tr g A is the mean curvature of X , defined with respect to the metric X*g where A is the 
second fundamental form, the mean curvature of X<p, and Cx* g '■= A 3 + \ A\ 2 . 

As mentioned, away from the central spheres, we need a scaled version in order to satisfy the 
hypotheses. We first prove a local estimate. 

Proposition A. 4. Let D C M be a disk of radius 1 in the X metric, centered at x. Let Y d £ : M 

R 3 be the immersion from Definition 4- 28. If v € C 2 '^(D, X ) satisfies 

ll»:^, D ,X)ll<^ 

then 

\\p- 2 (H v - H dX ) - C x v : C°'P(D,x)\\ < C( Cl )p(x)\\v : C 2 '? (D , X )\\ 2 . 

Proof. To apply Proposition A. 3, we rescale the target space by p{x). Thus, recalling item (8) 
from Lemma 4.44, the conditions (A.l) are satisfied and the hypothesis on the norm is satisfied for 
the dilated function p(x)v. Observe that H p ^ v = ~pj^H v and C(p( x )Y dX )*(fP(x)g) = p{ x Y' ''£(Y dX )* g - 
Thus, the proposition implies 

\\p~\x) (H v - H dX - C x « g v) : C°'"(AX)II < C( Cl )\\p(x)v : C 2 ^(D, X )\\ 2 . 

By the multiplicative property of the norms involved and again using item (8) of Lemma 4.44, we 
observe 

\\p- 2 (H v - H dX - C x * g v) : C°>P(D, X )\\ < C( Cl )p(x)\\v : C 2 '? (D , X )\\ 2 . 

□ 

To apply the quadratic part of the proposition in the main theorem, we need a global form. 
While better estimates can be shown to hold, we state the proposition in a way that is sufficient 
for our use. 

Proposition A. 5. Given a G (0, 1), let v G C 2 '@(M, R 3 ) such that |M| 2 ,/3,7 < r 1 ""/ 14 . Then 

\\p- 2 {H v - H dX ) - C x v\\ ,i3„ < T a - l \\v\\\ pa . 
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Proof. Notice that the definition of the global norm and of the function p implies that the hypothesis 
of Proposition A. 4 is satisfied for each x G M. The inequality from Proposition A. 4 implies that 
for any weighting /, 

\\p- 2 (H v -H d , ( -£ x . g v) : C°'P(D, X ,f)\\ < C(c 1 )f(x)p(x)\\v : C 2 ^(D, X , f)\ \ 2 - 

We now recall that the global norm implies that for Si \p] , Si [p, e, n] , we have / = 1 and for A\p, e, n] , 
f = e~ 7 -. Moreover, p is uniformly bounded above and below on each Si[p] and Si[p,e,n] for n 
even. Elsewhere, p(x) closely follows r(x) _1 (up to some small reparameterization) . Considering 
the various regions separately, we first observe that for D C Si [p] , 

C(ci)f(x)p(x)\\v : C 2 ^(D,x,f)\\ 2 < C\\v : C 2 ^(D, X )\\ 2 < r a ' l \\v\\l^ 

for any < a < 1 and r chosen small enough. The same inequality holds for Si\p, e,n] when n is 
even. Now consider that for n odd, p(x) < Ct^ 1 and / = 1. Thus for D C Si[p, e,n] with n odd, 
and < a < 1, 

C(ci)f(x)p(x)\\v : C 2 ^(D, X ,f)\\ 2 < °\\v : C 2 ^(D, X )\\ 2 



C 



< - T 2a -^ n \\v : C 2 'P(Si\p,e,n], X )\\ 2 < r^ 1 ] 



2 



Finally, we consider D C A[p, e,n']. Again, when n is even or odd, the problem changes. For n 
even, both / and p are decaying as functions of x and thus /(x)p(x) < Ct^ 1 for all x G A[p, e, ra'] 
for n even. But in this case 

C(c 1 )f(x)p(x)\\v : C?>P(D, X ,f)\\ 2 < °\\v : C 2 '@(D, x, e _7 -)|| 2 

< ^r 20 - 2 ^"- 1 )^ : (7 2 ^(Ab,e,n'], X ,e-^)|| 2 < t^IMI^. 

(Notice here we've been helped because n > 2.) If n is odd, observe that while / is decaying in x, 
/9 grows. But we note 



% V P r (e) / 

For < x < ^a, the control on the ratio p Td(e) (coming from (4.20)) and the fact that — 1 < w' < 

implies //? is decreasing as a function of x. Thus, f(x)p(x) < C(b) on all of A[p, e, n'] for n odd. It 
follows 



C{ci)f{x)p{x)\\v : C 2 ^(AX,/)|| 2 < C||u : C 2 ^(Ax,e~^) 



2 



< t"- 1 " 2 *- 1 ^ : C 2 ^(A[p,e,n'], X ,e-^)|| 2 < r^MI^. 

□ 

Appendix B. Dislocation at the Central Spheres 

For e G E(r) U i?(r), let /+[e] : M[e] n ([a, a + 3] x S 1 ) -> R and for e G £7(r) let /"[e] : 
M[e] n ([P[e] - (a + 3), P[e] - a] x S 1 ) ->• R such that on these regions 

(B.l) Ufe]- 1 o y d;C o D + = (Y + C\p + [e],e]) f+[e] 

(B.2) Ufe]- 1 o y d>c o £)- = (y + C[p-[e], e]) r[e] 

(B.3) u[e]- 1 oy diC o jD = (y + Cb + N,e]) /+[e] 

where D, are small perturbations of the identity map. Such functions exist, for r small, by 
Proposition A. 3. We prove important estimates for / + [e], e G E(T) and note that the same 
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estimates hold for the other two types of functions once we account for appropriate changes to the 
domain and notation. Notice that in the lemma below we use the convention of Section 6 and let 
C be a constant that may depend on b, f3, /?', 7, 7'. 

Lemma B.4. For f + [e] and e G E(T) as described above, we have the following: 

(1) f+[ e ] = on M[e] n ([a,a+l] x S 1 ). 

(2) ||/+[e] : C^(M[e] n ([a, a + 3] x S 1 ) < C||CII- 

(3) \\C x f+[e] - p- 2 H dislocation [e] : C°^(M[e] n ([a, a + 3] x S 1 ) , X )\ I < C||C|| 2 . 
(^J For /j as described in (5.16), 

3 

\\f + [e](x) +7r J2Q[P + [e], e]h(x) : C^(M[e] n ([a + 2,a + 3] x S 1 ) < C||C|| 2 . 
i=i 

Proof. First observe that on the domain of interest 

U[e]- 1 oy d)C = y edfle [r (e),r (i ( e ),/(e),C[p + [e],e],C[p-[e],e]]. 

Item (1) follows from the definition of the immersion and (2) follows from Proposition 3.28 once we 
note the domain of interest and the uniform equivalence of g, x on this domain. The local estimates 
of Proposition A. 3 and the C 2 ^ bound on f + [e] then imply 

\\C 90 f + [e]- p- 2 H dlslocation [e]oD + : C°'"(([a, a + 3] x S 1 ) n M[e], X )\ I < ^||C|| 2 . 

Here C go := A§2 + 2. We prove (3) by applying an appropriate modification Proposition 3.28. 
Finally, recalling the linear bound given in Proposition A. 3, and substituting Yo + £[p + [e], e] := Y[e] 
for X we observe 

||/ + [e]iV y[e] + C\p + [e], e] : C^(([a + 2, a + 3] x S 1 ) n M[e], X )|| < C||C|| 2 . 

Again using a modification of Proposition 3.28 and the fact that on ([a + 2, a + 3] x S 1 ) n M[e], 
Ufe]- 1 ^ = 7r(/i, h, /s), we prove (4). □ 

Definition B.5. For ease of notation, let 

S x [p] := M\p] U {j>=p+[e . ]} (M[ ej ] D [a, or] x S 1 ) 

U {j >= p - [e ,] } (M[ ej ] n [P[ ej ] - x,P[ ej ] - a] x S 1 ) . 

For each p € V(r), and {ej} = E p , let 

l\p]:S a+3 [p]^R 

such that 

' 0, ifx€M]p], 
f\p](x)=\ f + [e 3 ](x), ifp = p+[e J ],xeM[e j ]C\[a,a + ?>]x§\ 

k i/p = p-[e J -],x€M[e j ]n[P[e J -]-(a + 3),P[e J -]-a] x S 1 . 

The definitions of f\p] and ii dislocation immediately imply the following corollary. 

Corollary B.6. (1) f_[p] = on S a+1 [p\. 

(2) \\C x f[p] - p- 2 H dlslocati0 n : C°>^+ 3 [p], X )|| < CHCII 2 . 

(3) ifp = p + [e), \\M+irE-=i^lP + le},e}Mx) : C^(M[e]n{[a + 2, a + 3] x S 1 ) < C||C|| 2 
and 

(4) tfP = P~[e], 

3 

\\f[p] +vr^0[P + N,e]/ i (x) : C 1,/5 (M[e] n ([P[e] - (a + 3),P[e] - (a + 2)] x S 1 ) , x )|| < C||C|| 2 . 



i=l 
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